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ABSTRACT 

We use dipole and quadrupole statistics to test the large-scale isotropy of the first 1005 gamma-ray 
bursts observed by the Burst and Transient Source Experiment (BATSE). In addition to the entire 
sample of 1005 gamma-ray bursts, many subsets are examined. We use a variety of dipole and quadru- 
pole statistics to search for Galactic and other predicted anisotropies and for anisotropies in a 
coordinate-system independent manner. We find the gamma-ray burst locations to be consistent with 
isotropy, e.g., for the total sample the observed Galactic dipole moment <cos 0} differs from the value 
predicted for isotropy by 0.9 a and the observed Galactic quadrupole moment <sin 2 b — by 0.3 a. We 
estimate for various models the anisotropies that could have been detected. If one-half of the locations 
were within 86 of the Galactic center, or within 28° of the Galactic plane, the ensuing dipole or quadru- 
pole moment would have typically been detected at the 99% confidence level. We compare the observ- 
ations with the dipole and quadrupole moments of various Galactic models. Several Galactic gamma-ray 
bursts models have moments within 2 a of the observations; most of the Galactic models proposed to 
date are no longer in acceptable agreement with the data. Although a spherical dark matter halo dis- 
tribution could be consistent with the data, the required core radius is larger than the core radius of the 
dark matter halo used to explain the Galaxy’s rotation curve. Gamma-ray bursts are much more iso- 
tropic than any observed Galactic population, strongly favoring but not requiring an origin at cosmo- 
logical distances. 

Subject headings: cosmology: observations — gamma rays: bursts — methods: statistical 


1. INTRODUCTION 

The presence of a large-scale pattern in the locations of 
gamma-ray bursts (GRBs) would be a major clue to their 
origin. However, the gamma-ray burst locations determined 
by the Burst and Transient Source Experiment (BATSE) on 
the Compton Gamma Ray Observatory (CGRO) remain con- 
sistent with large-scale isotropy. This indicates either that 
GRBs are isotropic or that their anisotropy is too small to 
have yet been detected. The upper limits on the magnitudes 
of anisotropies strongly constrain Galactic models. 

One possible reason for not detecting an anisotropy is 
that we are observing gamma-ray bursts to a depth much 
less than the length scale of their distribution. This was 
thought to be the case prior to BATSE, when it was gener- 
ally believed that gamma-ray bursts originate from a disk 
population of neutron stars (Hurley 1986; Higdon & 
Lingenfelter 1990; Harding 1994). However, this possibility 
is excluded by BATSE’s concurrent observation that 
gamma-ray bursts are inhomogeneous, i.e., there is a defi- 
ciency of faint GRBs compared to a uniform density in 
Euclidean space. The inhomogeneity of the GRBs is demon- 
strated by <( F/K mllI > <0.5 and by the deviation of the 
log (N > P) — log P distribution from a -(3/2) slope 
power law for the fainter GRBs (Meegan et al. 1992b; 
Fishman et al. 1994; Meegan et al. 1994b; Horack & Emslie 


1994). The effect is not small: the number of GRBs with 
peak flux P = 0.8 photons s l cm 2 between 50 and 300 
keV on the 256 ms timescale is deficient by a factor of about 
6 relative to the extrapolation of the —(3/2) slope power- 
law characteristic of the homogeneous GRBs brighter than 
10 photons s 1 cm 2 (Pendleton et al. 1996). Assuming a 
Galactic origin, the observed inhomogeneity means that 
instead of viewing to less than a length scale, BATSE is 
viewing many length scales and sees past the “edge” of the 
distribution. 

Another possible explanation for not detecting an anisot- 
ropy, even though we are seeing past the edge of the dis- 
tribution, is that our offset from the center of the 
distribution is very small. A local version of this possibility 
is that gamma-ray bursts originate from a cometary cloud 
around the Sun (Bickert & Greiner 1993; Katz 1993; White 
1993; Luchkov 1994). In this hypothesis, the Earth-Sun dis- 
tance is too small relative to the diameter of the comet 
cloud to cause a detectable anisotropy and the edge of the 
comet distribution causes the inhomogeneity. A distant 
version of this possibility is that gamma-ray bursts orig- 
inate from a Galactic halo, in which case the length scale of 
the halo must be large compared to the distance from the 
Sun to the center of the Galaxy (Fishman 1979; Jennings 
1982; Shklovskii & Mitrofanov 1985; Atteia & Hurley 
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1986). Finally, gamma-ray bursts might be isotropic 
because they originate at cosmological distances (Prilutski 
& Usov 1975; Usov & Chibisov 1975; van den Bergh 1983; 
Paczynski 1986). In this hypothesis, the observed inhomo- 
geneity is due either to the non-Euclidean geometry of the 
universe, to source evolution, or to both. If GRBs come 
from z < 1, too few have been observed to detect the dipole 
moment due to our motion with respect to distant matter 
(Brainerd 1996). 

Arguably, the most crucial question about gamma-ray 
bursts is their distance. While discovering or placing strin- 
gent limits on large-scale anisotropies promises to disclose 
the distance scale, discovering gamma-ray burst repeaters 
would not, since there are cosmological models allowing or 
predicting repetition (McBreen & Metcalfe 1988; Dermer & 
Schlickeiser 1994; McBreen, Plunkett, & Metcalfe 1993; 
Brainerd 1994). The discovery of GRB repetition would, 
however, disprove the most popular model for a cosmo- 
logical origin, the mergers of compact stars. 

The BATSE team has described its findings that gamma- 
ray bursts are isotropic and inhomogeneous in previous 
papers and circulars (Meegan et al. 1991, 1992a, b, 1993a, 
1994b; Briggs et al. 1993a, b, 1994; Horack et al. 1993, 1994; 
Fishman et al. 1994). Many authors have emphasized that 
the isotropic and inhomogeneous GRB spatial distribution 
observed by BATSE is natural if gamma-ray bursts orig- 
inate at cosmological distances (Paczynski 1991a; Dermer 
1992; Mao & Paczynski 1992a; Piran 1992; Fenimore et al. 
1993; Wickramasinghe et al. 1993; Woods & Loeb 1994). 
Paczynski & Xu (1994) state that "A conservative conclu- 
sion is that the sources are at cosmological distances.” 
Other authors have emphasized the possibility of a Galactic 
origin with an anisotropic pattern too small to have been 
detected with BATSE, such as a very extended halo dis- 
tribution or a combination of source distributions (Brainerd 
1992; Eichler & Silk 1992; Hartmann 1992; Li & Dermer 
1992; Lingenfelter & Higdon 1992; Atteia & Dezalay 1993; 
Fabian & Podsiadlowski 1993; Liang & Li 1993; Smith & 
Lamb 1993; Hartmann et al. 1994; Higdon & Lingenfelter 
1994; Li, Duncan, & Thompson 1994; Lyne & Lorimer 
1994; Podsiadlowski, Rees, & Ruderman 1995). These 
authors have generally argued that GRB spectral and tem- 
poral properties make a Galactic origin likely. Quashnock 
& Lamb (1993a) identified a subset of the first BATSE 
catalog (IB; Fishman et al. 1994), which has a concentra- 
tion toward the Galactic plane and center, a pattern that 
they interpreted as demonstrating that all gamma-ray 
bursts originate in the spiral arms of the Milky Way. 

In this paper we test the isotropy of the locations of a 
larger sample, the first 1005 gamma-ray bursts detected by 
BATSE. Isotropy means that the probability of a GRB 
occurring in a given region of the sky is solely proportional 
to the solid angle of that region. While isotropy requires the 
absence of patterns on all angular scales, herein we empha- 
size the large angular scale properties of gamma-ray bursts 
by testing isotropy with dipole and quadrupole tests. We 
will analyze many subsets, including ones generated by the 
criteria of Quashnock & Lamb (1993a). We show that the 
locations are consistent with isotropy, derive upper limits 
on the parameters of several anisotropic models, and 
compare the observed moments with those of Galactic 
models. 

The most plausible anisotropy to which dipole and quad- 
rupole tests are insensitive is small-scale clustering, includ- 


ing GRB repetition, which averages to isotropy on large 
angular scales. Tests sensitive to such anisotropies include 
those based upon the nearest neighbor distribution and the 
two-point angular correlation function. Meegan et al. 
(1995), Hartmann et al. (1996), and Brainerd et al. (1995) 
apply these tests and others to the data of the second 
BATSE catalog (2B; Meegan et al. 1994a) and find no evi- 
dence for small-scale anisotropies in BATSE’s GRB loca- 
tions. While with the 2B catalog these authors do not 
confirm the evidence for GRB repetition found by Quash- 
nock & Lamb (1993b) and by Wang & Lingenfelter (1993, 
1995) in their analyses of the IB catalog, they do not com- 
pletely exclude burst repetition: the 99% confidence upper 
limit placed on the repeater fraction by Meegan et al. (1995), 
Hartmann et al. (1996), and Brainerd et al. (1995) is 20%. 

2. DIPOLE AND QUADRUPOLE TESTS AND DISTRIBUTIONS 

We test the large-scale isotropy of gamma-ray bursts by 
calculating various statistics that measure deviations from 
isotropy and comparing the observed values to the distribu- 
tion of values expected if GRBs were isotropic. If a statistic 
is found to have a value highly improbable assuming GRBs 
to be isotropic, then the null hypothesis of isotropy is sta- 
tistically contradicted. Several different statistics are used 
because each is most sensitive to finding certain patterns of 
anisotropies. 

Large angular scale anisotropies are best searched for 
with dipole and quadrupole statistics. Various dipole and 
quadrupole statistics have been introduced to the field of 
gamma-ray burst studies by Hartmann & Epstein (1989), 
Paczynski (1990), and Briggs (1993) and are reviewed by 
Briggs (1993). Dipole statistics are sensitive to a concentra- 
tion of GRBs toward one direction on the sky, while quad- 
rupole statistics are sensitive to concentrations in a plane or 
toward two opposite poles. Any plausible large-scale con- 
centration should be revealed by its dipole or quadrupole 
moments before higher moments become significant. 

Dipole and quadrupole statistics based upon a particular 
coordinate system are used since they will be most sensitive 
to anisotropies in that coordinate system. Galactic-based 
tests are used because they are most sensitive to Galactic 
anisotropies, which would result from a Galactic origin for 
GRBs. Equatorial-based statistics are used because these 
statistics are most sensitive to the artificial anisotropy 
induced by BATSE’s nonuniform sky exposure. We calcu- 
late the dipole moment toward the Sun and the quadrupole 
moment in the ecliptic plane to test for a heliocentric origin. 
A final pair of dipole and quadrupole statistics are con- 
structed in a coordinate-system independent manner in 
order to search for anisotropies in unexpected directions in 
a model-independent manner (Briggs 1993). These statistics 
are characterized in Table 1. 

Each statistic is used as a one-tailed test : the sign of the 
deviation that is significant is indicated in the last column of 
Table 1. We use the coordinate-system based tests as one- 
tailed tests because the models tested, a Galactic origin, a 
solar system origin, and BATSE’s sky exposure, predict the 
sign of the expected deviation. We use the coordinate- 
system independent tests iK and ^ as one-tailed tests 
because no random anisotropic distribution can cause small 
values of these statistics. 

We use two probability distributions on the sphere in 
Monte Carlo simulations designed to test the methods. A 
standardized form of the Fisher distribution has a dipole 
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TABLE 1 

The Dipole and Quadrupole Statistics 


Statistic* 

Coordinate 

System 

Moment 

Tested 

Asymptotic 

Distribution 

Mean /x 
for Isotropy 

a 

Sensitive to Predicted 
Concentration at 

Signature of 
Predicted 
Concentration 

<cos 0) 

Galactic 

Dipole 

Gaussian 

0 

Jvin b 

Galactic center 


<sin 2 b — 

Galactic 

Quadrupole 

Gaussian 

0 

V4/45N. 

Galactic plane 


Rayleigh-Watson # 

Independent 

Dipole 

xi 

3 


Any point 


Bingham M 

Independent 

Quadrupole 

, xi . 

5 


Any plane or 2 poles 


<sin <5) 

Equatorial 

Dipole 

Gaussian 

0 

Vi/ 3/v i> 

Earth’s N pole 


<sin 2 S — 

Equatorial 

Quadrupole 

Gaussian 

0 

•J 4/45N B 

Earth’s poles 


<cos 0) 

Heliocentric 

Dipole 

Gaussian 

0 

V1/3JV, 

Sun 

>V 

<sm-’/f J> 

Ecliptic 

Quadrupole 

Gaussian 

0 

V 4 /45«. 

Ecliptic plane 

<n 


a 0 is the angle between a burst and the Galactic center, b is Galactic latitude, 6 is declination, 4> is the angle between a burst and the Sun, and ft is ecliptic 
latitude. 


moment toward the Galactic center — the probability 
density per solid angle is proportional to exp (tc cos 0\ 
where 0 is the angle between a point and the Galactic center 
(Fisher, Lewis, & Embleton 1987). Increasing values of the 
concentration parameter k result in a greater concentration 
toward the Galactic center (Fig. la). A standardized form of 
the Watson distribution has a Galactic-based quadrupole 
moment — the probability density per solid angle is pro- 
portional to exp (k sin 2 b\ where b is Galactic latitude 
(Fisher et al. 1987). Values of the concentration parameter 
k < 0 result in a concentration about the Galactic plane, 
with more negative values of k resulting in greater concen- 
tration (Figs, \b-\d). Values of k > 0 result in a concentra- 


tion at the Galactic poles. Fisher et al. (1987) and Best & 
Fisher (1986) give algorithms for simulating locations from 
these distributions. 1 

Figure 2 characterizes these distributions, showing, as a 
function of k , the moments of the distributions. Also shown 
for the Fisher distribution is the angle 0 1/2 such that the 
region within d l/2 of the Galactic center contains one-half of 

1 Both references have typographic errors in the algorithm for the 
Watson distribution. The algorithm of Fisher et al. (1987) should have 
“0 - arccos S” instead of “0 = 5”. Also note that this algorithm gener- 
ates locations only for one hemisphere— locations must be moved to the 
other hemisphere with probability one-half. 



Fig. 1 a 


Fig. 1 h 


f 90 


Fig. lc 


1 HO 



Fig. Id 


Fig. 1. Simulations of 1000 locations from several anisotropic distributions shown on Hammer- Aitoflf projections. The simulations are typical in 
the sense of having moments within 20% of the expected value, (a) Fisher distribution with k = 0.13. This distribution has one-half of the probability 
within t) v 2 = 86 of the Galactic center and <cos 0 > = 0.043. ( b ) Watson distribution with k = -0.28. This distribution has one-half of the probability within 
b v2 = 28 "of the Galactic plane and <sin 2 b - \) = —0.024. (c) Watson distribution with k = — l. This distribution has one-half of the probability within 
b v2 = 23 of the Galactic plane and <sin 2 b - = -0.080. (d) Watson distribution with k = - 10. This distribution has one-half of the probability within 

b l/2 = 8 7 of the Galactic plane and <sin 2 b - = -0.283. 
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Fig. 2. — Properties of the Fisher and Watson distributions vs. concen- 
tration parameter k. (a) Fisher distribution: 0 l/2 such that P(0 < 0 17 ) = 
0.5 {solid curve ) and dipole moment (cos 0 ) ( dashed curve); ( b ) Watson 
distribution: b lf2 such that P{\b\ < b ll2 ) = 0.5 {solid curves) and quadru- 
pole moment (sin 2 b — -j) {dashed curves ). The upper curves are for the 
bipolar case, k > 0, while the lower curves are for the plane-concentration 
case, k; < 0. 

the probability. For the Watson distribution, Figure 2b 
shows the angle b lf2 such that the region within b Xj2 of the 
Galactic plane contains one-half of the probability. 

3. ANALYSIS 

3.1. Location Determination 

BATSE consists of eight detector modules located at the 
corners of the Compton Gamma Ray Observatory (CGRO). 
Each module contains a Large Area Detector (LAD) and a 
Spectroscopy Detector (SD). The results presented here are 
based upon the LADs, which are 50.8 cm diameter by 1.27 
cm thick Nal(Tl) scintillators. The instrument triggers and 
records extensive information whenever there is a 5.5 o 
increase about background in two or more detectors. 
Further descriptions of the instrument and the various data 
types transmitted are given by Horack (1991) and by 
Fishman et al. (1994). 

BATSE determines GRB locations by comparing the 
rates in the eight LADs (Brock et al. 1992b; Fishman et al. 
1994). The locations are determined by the program LOC- 
BURST by fitting the signal in four detectors and mini- 
mizing x 2 - The model includes the direct response of the 
detectors and scattering in the spacecraft and from Earth. 

The locations are not precise because of both statistical 
and systematic limitations. The statistical limitations are 
due to Poisson fluctuations in the counts observed with the 
detectors — these errors <r stat are determined by LOC- 
BURST from the derivatives of / 2 , which are obtained from 
the observed counts and the detector model. These errors 
are well understood and are Gaussian in the count regime 
of BATSE. 


Normally there are several data types to choose from in 
order to obtain the best possible location, but due to the 
failure of the CGRO tape recorders, for 121 GRBs of our 
total sample of 1005 GRBs only the MAXBC data type is 
available, which consists of the count rates in each detector 
for the peak 1 s of emission. Because the background is 
determined differently for MAXBC data than for other data 
types, until very recently LOCBURST was unable to esti- 
mate statistical errors for M AXBC-determined locations. 

The statistical errors range from a fraction of a degree for 
the brightest GRBs to typically 13° for GRBs at the trigger 
threshold, as summarized in Table 2. Because the (re- 
values were unavailable for MAXBC-located GRBs, these 
GRBs are not included in Table 2. For bursts in which the 
MAXBC data type must be used to determine the location, 
somewhat larger statistical errors are expected: (1) for 
bursts longer than 1 s, because of the neglect of source 
counts outside the peak 1 s interval, and (2) for bursts 
shorter than 1 s, because of the reduced SNR caused by 
including background-only times in the 1 s long source 
interval. After most of the analysis of this paper was accom- 
plished, LOCBURST was improved to correctly estimate 
the statistical errors of MAXBC-determined locations. 
Based upon the currently available sample of 24 repro- 
cessed MAXBC locations, the statistical error distribution 
of MAXBC-determined locations is only moderately worse 
than that of locations determined with other data types: 14 
of the 24 GRBs have a stat < 3?7 and only three have <r slat > 
13?8, with the largest value being 24°. 

The accuracy of BATSE’s GRB locations have been veri- 
fied in several ways. For bright GRBs accurately deter- 
mined locations are available from the Interplanetary 
Network (IPN; Hurley 1993). Comparison of these loca- 
tions with the locations determined from the LAD rates 
show that the systematic error of the BATSE locations is 
currently typically 4° (Fishman et al. 1994; Hurley et al. 
1994). These systematic errors are due to calibration uncer- 
tainties, imperfect background subtraction, and other 
approximations in the data analysis. We estimate the 1 o 
location error for a GRB as the root mean square sum of its 
estimated statistical error and a 4° systematic error. Since 
the 4° error is a systematic error, its distribution need not be 
Gaussian, and we find from comparison of LAD- 
determined locations with IPN-determined locations that 
the total error is approximately Gaussian, with somewhat 
larger “tails” (Fishman et al. 1994). The excess tails of the 
distribution are due to the systematic error component and 
are unimportant for bursts for which the statistical error is 
dominant. 

In addition to their somewhat larger statistical errors 
<r stat , MAXBC-located GRBs should have somewhat larger 
systematic errors o sys because the background is determined 


TABLE 2 

Angular Errors: Data Set 6 and Simulations 



2B LESS MAXBC 

Simulations 

Fraction 



Fraction 

I 

50% 

0? 1 —3.7 

4.0-5.4 

45% 

15: 

25 

3. 7-7. 4 

5.4-8. 4 

22.5 

20. 

18.8 

7.4-13.8 

8.4-14.4 

16.875 

30. 

6.2 

13.9-29.9 

14.5-30.1 

15.625 

60. 
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on-board the spacecraft by a simpler algorithm than that 
used by LOCBURST on the ground. Koshut et al. (1994) 
analyzed the MAXBC-determined locations of 33 events 
with locations determined by the IPN or known from their 
solar origin and found that the systematic error was < 7°. 

The locations have also been verified by examining events 
from sources of known location: solar flares, Cygnus X-l 
fluctuations, and triggers from SGRs (Fishman et al. 1994; 
Meegan et al. 1993b; Kouveliotou et al. 1993a, 1994). These 
verifications are important since they are not limited to the 
bright events required for IPN -locations: solar flares can be 
of any intensity, and Cygnus X-l fluctuations and SGRs 
events are near threshold. Fluctuations from Cygnus X-l 
are more difficult to locate than GRBs of comparable inten- 
sity because they are always short and are superposed on a 
background of subtrigger fluctuations. The distribution of 
separations between the location of Cygnus X-l and the 
LOCBURST-determined locations of the fluctuations has a 
standard deviation of 13° (Meegan et al. 1993b), the typical 
value of cr stat for events near threshold. For these weak trig- 
gers, the systematic error er sys = 4° is unimportant. 

The program LOCBURST has undergone steady 
improvement: the version used to determine the post-2B 
locations used herein is slightly better than the version used 
for the post- IB locations of the 2B catalog, which is dis- 
tinctly better than the version used for IB locations. After 
the calculation of the locations used herein, major improve- 
ments have been made to LOCBURST, which reduce the 
systematic error to <2°. As we demonstrate below, dipole 
and quadrupole moments are relatively insensitive to loca- 
tions errors so that all of the locations are useful for this 
analysis. The location improvements are important for tests 
of clustering or repetition, since the statistical signal rapidly 
strengthens with improved locations (Brainerd et al. 1995; 
Hartmann et al. 1996). At the time of this writing 95 of the 
locations used herein, all post-2B, have been redetermined 
with the newest version of LOCBURST. Comparison of 
these revised locations with the locations used herein show 
that 50% of the locations have changed by less than 4?8 and 
only 10% have changed by more than 14?0. The largest shift 
is 28°. In terms of the old total error estimates a tol = 
(4 2 + <r slat ) 1/2 , 50% have moved by less than 0.9a tot , and 
90% have moved by less than 2.7a tol . The worst shift in 
terms of a tot is an 18° change that is 4.6a tot . 

The anisotropic response of BATSE's detectors, which 
enables determination of GRB locations, implies an inten- 
sity threshold for detecting a GRB, which varies across the 
sky (Brock et al. 1992a). This intensity-dependent, nonuni- 
form sky response is best described in CGRO coordinates. 
Because of the many orientations of CGRO for the various 
observation periods, this effect averages away and need not 
be further considered. 

Because CGRO is in a low-Earth orbit, about one-third of 
the sky is blocked by the Earth. This causes the equatorial 
regions to be observed with about 20% less exposure than 
the polar regions. Additionally, because burst triggering is 
disabled during passages through the South Atlantic 
Anomaly (SAA) and several other regions of enhanced 
charged particle background, the south pole receives about 
15% less exposure than the north pole. Both of these effects 
are intensity independent and are best described in equato- 
rial coordinates: right ascension, a, and declination, 3. Over 
the spacecraft orbital precession period of about 50 days, 
the right ascension dependence of the exposure averages 


away, leaving the sky exposure a function only of decli- 
nation: T(3 ). The sky exposure T(<5) has been measured 
(Fishman et al. 1994), allowing corrections for this effect to 
be made. 

3.2. Analysis Methods 

The values of the statistics are calculated from the loca- 
tions; consequently, if many samples are 44 drawn” from a 
population, the statistics will have different values for each 
of the samples. The fluctuations in a finite sample limit the 
ability of even a perfect instrument to detect a small anisot- 
ropy. To demonstrate this effect, Monte Carlo simulations 
of 250 isotropic locations were made. Figure 3 shows histo- 
grams of the values of <cos 0> and <sin 2 3 — obtained 
from the simulations overlaid with the theoretically 
expected asymptotic distributions (solid curves). As a char- 
acterization of the expected fluctuations, Table 1 lists the 
expected standard deviations a of the six statistics, derived 
from the asymptotic distributions. The asymptotic distribu- 
tions may not be sufficiently accurate for small numbers of 
gamma-ray bursts N B and extreme values of the statistics, 
which is especially likely for <sin 2 b — (Briggs 1993). 
When necessary we use the actual distributions of the sta- 
tistics as determined by Monte Carlo simulations of isot- 
ropy. 

An artificial anisotropy should exist in the observed 
sample of GRB locations because BATSE’s sensitivity is not 



-.1 -.05 0 .05 .1 

(cos 6) 


CD 



-.05 „ 0 05 

(sin <5 - 1/3) 


Fig. 3. — Simulations of 250 isotropic locations with and without inclu- 
sion of the effects of BATSE's sky exposure. The histograms overlaid with 
solid curves are simulations of the ideal case. The solid curves are the 
expected asymptotic distribution of the statistics. The histograms overlaid 
with dashed curves are the result of including BATSE’s sky exposure map 
in the simulations. The dashed curves are Gaussians of the same means 
and standard deviations as the sky-exposure biased histograms, (a) Galac- 
tic dipole moment <cos 0 >. The sky-exposure induced shift is minor 
because the artificial dipole is moderate and only part of it projects toward 
the Galactic center, (b) Ecliptic quadrupole moment <sin ;! <5 — -j). This 
statistic detects the shadowing of the GRBs by the Earth; it has a relatively 
large sky-exposure induced shift. 
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uniform across the sky. The resulting bias expected in the 
statistics has been found by Monte Carlo simulations — 
isotropic locations are created using a random number gen- 
erator and included in the simulated sample with a 
probability obtained from the sky exposure map T(8). 
Locations are generated until N B have been accepted and 
then the values of the statistics are calculated. This effect is 
illustrated in Figure 3 by the second histograms with their 
overlying dashed curves. The histograms are still Gaussian, 
but the artificial anisotropy of the sky exposure has caused 
the expected mean of the observation to shift from zero. In 
addition, there are very small changes in the widths of the 
distributions. The shift in the mean is large for <sin 2 3 — y> 
because this equatorial-based statistic is sensitive to the 
time-averaged shadowing of the celestial equator by the 
Earth. The shift in the mean of <cos 0 ) is small because the 
dipole induced by the equatorial north-south asymmetry is 
a smaller effect, relative to the standard deviation o of the 
fluctuations, than the quadrupole moment induced by the 
Earth blockage and because only part of the equatorial 
dipole is projected onto the direction toward the Galactic 
center. The fraction of the artificial dipole moment project- 
ed equals the cosine of the angle between the equatorial 
plane and the Galactic plane, which is cos (60°) = 

Table 3 lists the means and standard deviations of the 
eight statistics expected for isotropy modified by BATSE’s 
actual sky exposure. These sky-exposure-corrected mean 
values of the statistics have been determined using BATSE’s 
sky exposure accumulated for the time period of the first 
BATSE catalog (Fishman et al. 1994), the 320 days from 
1990 April 21 to 1992 March 5. Work is in progress to 
measure the sky exposure after this period, a task made 
difficult by the failure of the CGRO tape recorders. The 
expected values have also been calculated for the sky expo- 
sure map for the first 165.5 days (Brock et al. 1992a) — the 
differences are insignificant. We expect the changes in 
BATSE’s sky exposure with time to remain small, so the 
extrapolation of the current sky exposure map to N B = 
2000 gamma-ray bursts should be a good approximation. 

These two effects, the finite sample fluctuations and the 
shifts in the expected means due to BATSE’s sky exposure, 
are illustrated in Figure 4. The dashed lines show the results 
expected ignoring BATSE’s sky exposure, while the solid 
lines indicate the results expected based upon BATSE’s sky 
exposure. The bold lines show the expected means of the 

TABLE 3 

Expected Values of the Statistics: Corrected for BATSE’s 
Sky Exposure 


Statistic Mean a,b 


<COS0> —0.013 0.99^1/3 N B 

<sin 2 b - -0.005 o.99 v /4/45N B 

IT 3 4* 0.0020N* 

£ 5 + 0.0077V fl 

<sin<5> 0.026 1.04yi/3N B 

<sin 2 S - 0 0.026 1.03^4/45 N~ B 

<cos 4>} 0.010 cos [2n(T - 9160)/365.25] 0.98^1/3/% 

(sin 2 /? - j) 0.020 1.02^/4/45 N B 


* These means are valid for time intervals long enough to average over 
the CGRO precession period of about 50 days. 

b T is the time in truncated Julian day number (TJD), which is the 
Julian day number minus 2440000.5 
c See Fig. 4. 


statistics, while the normal-weight lines show the + 1 a 
envelopes of the distributions of the statistics. The widths of 
the distributions (indicated by the outer lines) are due to the 
finite sample fluctuations, while the shifts of the means 
(indicated by the separation between the dashed, bold lines 
and the solid, bold lines) are due to the artificial anisotropy 
induced by BATSE’s sky exposure. 

The lines on Figure 4 allow visual evaluation of the effects 
of BATSE’s nonuniform sky exposure. We consider the 
effect of the nonuniform sky exposure to be major when it 
causes the corrected mean of a statistic to differ from the 
ideal mean by > 1 <r, i.e., when the bold, solid line is outside 
of the outer dashed lines. For the Galactic-based statistics, 
<cos 0) and <sin 2 b — |), this only happens for N B > 2000. 
These statistics are relatively insensitive to the nonuniform 
sky exposure because the Earth’s equator is highly inclined 
to the Galactic equator. The coordinate-system indepen- 
dent quadrupole statistic M is not importantly affected until 
N b > 400, while the coordinate-system independent dipole 
statistic iC is not importantly affected by the nonuniform 
sky exposure until N B > 1200. Because the nonuniform sky 
exposure originates in the equatorial coordinate system, the 
equatorial-based statistics are quite sensitive to it: the 
effects are important for <sin 3 > for N B > 500 locations and 
for <sin 2 3 — for N B > 130 locations. The dipole statistic 
<sin 3 ) detects the ~15% increased exposure of the north 
pole relative to the south pole, while the quadrupole sta- 
tistic <(sin 2 3 — detects the ^20% reduced exposure of 
the equator relative to the poles. 

As described above, the locations of the GRBs are uncer- 
tain, which propagates into uncertainties in the values of the 
statistics. Because the location errors are small compared to 
the scale of large-scale anisotropies, the corresponding 
errors in the values of the statistics are negligible. Monte 
Carlo simulations were conducted to demonstrate this. The 
histograms on the left of each panel of Figure 5 are based 
upon 5 x 10 5 simulations of 1000 isotropic locations. The 
curves overlaying these histograms are the asymptotic dis- 
tributions of the statistics. The vertical lines on Figure 5 
show the critical values of the statistics corresponding to 
confidence levels against the null hypothesis of isotropy of 
99%, 99.9%, and 99.99%. These critical values are defined 
by the requirement that the chance probabilities under the 
assumption of isotropy of more extreme values of the sta- 
tistics are 1%, 0.1%, and 0.01%. 

Simulations were made of two distributions with moder- 
ately detectable anisotropies. The results of the simulations 
made assuming perfect determinations of the locations are 
shown in the bold, rightmost histograms on Figure 5. The 
fractions of these histograms to the right of the vertical lines 
are the fractions of the simulations in which the aniso- 
tropies are detected at the confidence levels of the lines. The 
anisotropic simulations were repeated, this time moving 
each simulated location an angle E in a random direction, 
obtaining the middle histograms of Figure 5. 

The angles E and the fractions of the simulated bursts 
moved by each particular value of E are listed in the last 
two columns of Table 2. The simulated smearing exagger- 
ates with a very substantial safety margin the actual loca- 
tion uncertainties: (1) the values E are based upon the worst 
<7 tot in each row, (2) the values E are at least twice this worst 
location error, (3) the simulated smearing locations are each 
moved by an angle E, not by an angle selected from a 
Gaussian distribution of this width, and (4) to allow for 
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Fig. 4.— These panels show the mean values and standard deviations expected for isotropy of several dipole and quadrupole statistics for N - 20 to 2000. 
Also depicted are selected observed values of the statistics. The dashed lines indicate the expected values of the statistics for an ideal observation, while the 
solid lines indicate the values expected based upon BATSE’s sky exposure. The bold lines indicate the means of the statistics, while the normal-weight lines 
show the ± 1 a envelopes of the distributions of the statistics. Observations {dots) are shown for the first 40, 100, 262 ( = IB revised), 585 (=2B) and 1005 
gamma-ray bursts observed by BATSE. These are not independent sets. The observed values should be compared to the solid lines. When the observed 
values are within the normal-weight, solid lines they deviate by less than 1 o from the values expected for isotropy, (a) Galactic dipole statistic, <cos O ' >, ( b ) 
Galactic quadrupole statistic, <sin 2 b — (c) coordinate-system independent Rayleigh-Watson dipole statistic W , (d) coordinate-system independent 

Bingham quadrupole statistic if, ( e ) equatorial dipole statistic, <sin <5>, (/) equatorial quadrupole statistic, <sin 2 £>. 


outliers additional simulated locations are placed into the the Galactic center and has an expected dipole moment of 
category with the largest smearing angle. <cos 6 > = 0.06649. The middle histograms, which show the 

Figures 5a and 5b show simulations of the Fisher dis- values of the statistics when the locations are smeared, are 

tribution with k = 0.2, a slightly more concentrated dis- only slightly shifted from the bold histograms on the right, 

tribution than the example shown in Figure la. This which show the results for perfect locations. Because the 

distribution places one-half of the probability within 84?3 of scale of the anisotropy is of order 90°, the smearing by tens 
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l^l 5 ^' mUlatl ° nS W1 u h and Wit J h0Ul effects of exaggerated locat,on errors - The histograms on the left, overlaid with curves, are of 1000 isotropic 
qq qo'v fia CUrV | eS are . the expected asymptotic distributions The three vertical lines define the critical values corresponding to the 99%, 99.9%, and 
. «» confidence levels, i.e values of the statistics to the right of these lines are less probable under the isotropic hypothesis than 1 %, 0.1%, and 0 01 % The 

I™ r 6 nghl ’ Wh j C L arC b °' d ’ ar , e , f ° r slmulated anisotropic locations without any location errors. The middle histograms are for simulated 
the F Ih P H ^ K°? S smear . ed b i " 1 ° vlng f aL ; h location m a random direction by angles Z as specified in Table 2. For {<:) and (b) the anisotropic distribution is 
n u/ W h k a- °7’ Whl C f ° r ,C) and (d) , the dlstnbu,lon is mixed: 90% isotropic, 10% from the Watson distribution with k = - 10 (a) 
Bingham quadrupolsUbll ' P Sta “ S " C # 1 W <Sm ’ * _ 5> ’ n °' e ‘ hat ‘ he X ' aX ' S ' S reVCrSed SO ,hat ,he anisotropic histograms are to the right; (d) 


of degrees is a small effect. To test the effect of the location 
errors for a quadrupole distribution, simulations of the 
Watson distribution for k = — 10 (Fig. Id) were made — this 
distribution has one-half of the probability within | b \ < 8?7 
and an expected quadrupole moment of <sin 2 b — 3) = 
— 0.2833. This distribution is extremely detectable, so its 


detectability was reduced by diluting it with isotropic loca- 
tions, creating anisotropic locations with probability 10%. 
In this example (Fig. 5c and 5d ), the middle histograms 
based upon the smeared locations are quite noticeably 
shifted from the bold histograms based upon the perfect 
locations, but the smearing has not caused the middle histo- 
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grams to be similar to the histograms for isotropic loca- 
tions. 

The sensitivity or power of a statistic is the probability of 
its finding an actual anisotropy (Eadie et al. 1971; Martin 
1971). In our case, the powers of the statistics for finding the 
example anisotropies at the confidence levels of 99%, 
99.9%, and 99.99% is given by the fraction of simulations 
with statistics more extreme than the corresponding critical 
values, i.e., to the right of the vertical lines in Figure 5. The 
powers of the statistics are listed in Tables 4A-4D. The 
decrease in the fraction of the anisotropic simulations, 
smeared versus unsmeared, that are past the critical values 
represents the loss in sensitivity or power to the anisotropy 
caused by the smearing. Table 4A shows that the loss of 
power for the Fisher k = 0.2 model is minor. Similar results 
were obtained by Horack et al. (1993) via simulations of a 
halo distribution. The model with 90% isotropic locations 
and 10% locations from a narrow disk model was intended 
to be an extreme case — half of the disk probability is con- 


tained within | b | < 8?7 versus the exaggerated smearing 
angles E listed in Table 2. Table 4B shows that while the 
smearing substantially reduces the sensitivity to this model 
anisotropy, the ability to detect the narrow disk component 
is not eliminated. 

The above simulations were each of 1000 locations. We 
now discuss how our ability to find anisotropies with 
BATSE depends on the number of gamma-ray bursts N B 
observed. The finite-sample fluctuations of the coordinate- 
system based statistics, such as <cos 0> and <sin 2 b - j>, 
asymptotically have a l/(N B ) 112 dependence (Tables 1 and 3) 
because of the Central Limit Theorem (Paczynski 1990; 
Briggs 1993). Monte Carlo simulations show that the 
asymptotic distributions are sufficiently accurate for the 
values of N B under consideration (Briggs 1993). The 
coordinate-system independent Rayleigh- Watson H r and 
Bingham & statistics are proportional to the product of N B 
and the moment per GRB squared (Bingham 1974; Watson 
1983; Briggs 1993), and their asymptotic distributions are 


TABLE 4A 

5 X 10 5 Simulations of 10' Locations from Fisher Distribution with k = 0.2 




<cos 0 ) 



IT 





Simulations Past Confidence Levels 


Simulations Past Confidence Levels 

Locations 

Mean 

99% 

99.9% 99.99% 

Mean 

99% 

99.9% 

99.99% 

Perfect locations 

Smeared locations 

Factor change 

0.06647 ± 0.00002 
0.05786 + 0.00002 
0.87 

90.7% 

80.1 

71.0% 46.9% 

53.2 29.1 

16.24 + 0.01 
13.04 + 0.01 
0.80 

71.4% 

54.5 

44.8% 

27.9 

23.8% 

12.2 




TABLE 4B 






5 x 10 5 Simulations of 10 3 Locations, 10% from Watson Distribution with k = - 10.0 





<$in 2 b — 

i> 


3 





Simulations Past Confidence Levels 


Simulations Past Confidence Levels 

Locations 

Mean 

99% 

99.9% 99.99% 

Mean 

99% 

99.9% 

99.99% 

Perfect locations 

Smeared locations 

Factor change 

— 0.02834 ± 0.00001 
-0.01915 ± 0.00001 
0.68 

75.4% 

38.4 

46.7% 23.7% 

14.2 4.3 

14.02 ± 0.01 
9.12 + 0.01 
0.65 

38.7% 

12.5 

16.3% 

3.1 

5.9% 

0.7 




TABLE 4C 






5 x 10 4 Simulations of 10 4 Locations from Fisher Distribution with k - 

0.063095 





<cos oy 



IT 





Simulations Past Confidence Levels 


Simulations Past Confidence Levels 

Locations 

Mean 

99% 

99.9% 99.99% 

Mean 

99% 

99.9% 

99.99% 

Perfect locations 

Smeared locations 

Factor change 

0.02104 ± 0.00003 
0.01830 + 0.00003 
0.86 

90.7% 

80.0 

71.1% 47.1% 

53.4 29.4 

16.28 + 0.03 
13.05 ± 0.03 
0.80 

71.5% 

54.5 

45.1% 

28.1 

24.0% 

12.2 




TABLE 4D 






5 x 10 4 Simulations of 10 4 Locations, 3.162% from Watson Distribution with k = - 10.0 





<sin 2 b - 

i> 


3 





Simulations Past Confidence Levels 


Simulations Past Confidence Levels 

Locations 

Mean 

99% 

99.9% 99.99% 

Mean 

99% 

99.9% 

99.99% 

Perfect locations 

Smeared locations — 
Factor change 

-0.00896 + 0.00001 
-0.00607 + 0.00001 
0.68 

75.0% 

38.6 

46.7% 23.7% 

14.7 4.6 

14.01 ± 0.03 
9.14 + 0.02 
0.65 

38.4% 

12.7 

16.2% 

3.2 

5.9% 

0.7 
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also sufficiently accurate (Briggs 1993). Thus, for both types 
of statistics the finite sample fluctuations cause the magni- 
tude of detectable anisotropic moments to scale as 
*/(N b ) . The effects of the location uncertainties scale in 

the same manner. It does not matter where a variance 
comes from, only that it is finite, in which case the Law of 
Large Numbers or the Central Limit Theorem tells us that 
the mean of the measurements will be determined with 
uncertainty scaling as l/(N B ) 1/2 for large N B (Eadie et al 
1971; Martin 1971; Lyons 1986). This is the well-known 
result that the knowledge of the mean of a series of numbers 
with errors improves as l/(iV) 1/2 . 

As a demonstration that BATSE’s ability to find aniso- 
tropies improves as 1/(N B ) 1 ' 2 , 5 x 10 4 simulations of 10 4 
locations were made (Tables 4C and 4D). The simulations 
in Tables 4C and 4D have 10 times the number of locations 
as those of Tables 4A and 4B, but the moments have been 
reduced by a factor of (10) 1/2 . Allowing for the fluctuations 
expected in 5 x 10 4 simulations, the anisotropies are 
detected in the same fraction of simulations in Tables 4C 
and 4D as in Tables 4A and 4B, confirming the l/(N„) 112 
improvement in sensitivity. 

Smearing by angles Z as listed in Table 2 reduces the 
dipole moment <cos 6 > of dipole distributions by a factor of 
0.86 and the quadrupole moment <sin 2 b - of quadru- 
P°le distributions by a factor of 0.68. We found this to be 
true not only of the distributions listed in Tables 4A-4D, 
but also for all the other distributions we simulated. These 
reduction factors greatly overestimate the effects of 
BATSE’s location uncertainties since the angles Z greatly 
exaggerate the location uncertainties. 

Another possible error in the values of the statistics is due 
to the possibility of trigger misclassifications. We estimate 
both the number of GRBs misclassified as non-GRBs and 
the number of non-GRBs misclassified as GRBs to be about 
1% (Mallozzi et al. 1993; Fishman et al. 1994). The esti- 
mated misclassifications are roughly equally split between 
Cygnus X-l fluctuations and solar flares, so the misclassifi- 
cations are not all at the same location. We also estimate 
that at most one soft gamma-ray repeater (SGR) has been 
misidentified as a GRB and vice versa (Kouveliotou 1994). 
Even if all nine SGR triggers are added to the set of 1005 
gamma-ray GRBs that we analyze, the set remains consis- 
tent with isotropy. In an extreme example, if a false popu- 
lation of 2% of the GRBs is placed at the Galactic center, a 
dipole moment of <cos O') = 0.02 is produced, which would 
be detectable only at the 1.1 a level in a sample of 1000 
locations. 

In summary, the dominant limitation for BATSE in dis- 
covering a large-scale anisotropy is the fluctuations 
resulting from the finite number of GRBs observed, a limi- 
tation that no instrument can avoid. The effects of the loca- 


+90 



tl0nS TI 0f the firSt 1005 S amma -™y bursts observed by 
BATSE. The map is a Hammer-Aitoff projection in Galactic coordinates. 


tion uncertainties are negligible and will remain so since the 
effects of both the finite-sample fluctuations and the loca- 
tion uncertainties scale as l/(JV fl )‘' 2 The estimated rate of 
tngger misclassifications is too low to significantly effect the 
results. The bias induced by BATSE’s nonuniform sky 
exposure can be measured and corrected for when compar- 
ing model predictions with the observations. 


4. RESULTS 

The first 1005 gamma-ray bursts observed by BATSE are 
shown in Figure 6. These GRBs occurred between the 
enabling of burst triggering on 1991 April 21 and the end of 
CGRO Viewing Period 222 on 1994 May 31. The post-2B 
catalog locations are preliminary, but as discussed in § 3 1 
the revised locations will be close to the current locations 
and are adequate for the study of the large-angular scale 
properties of gamma-ray bursts. The 585 locations of the 
second BATSE catalog (2B) are publicly available (Meeean 
et al. 1994a). The first 262 triggers of the 2B catalog are 
referred to as the revised IB catalog— these differ from the 
original IB catalog (Fishman et al. 1994) by the revised 
classifications of several triggers. 

The values of the six statistics for the first 40 100 262 
(=1B revised), 585 (=2B) and 1005 gamma-ray bursts 
observed by BATSE are shown on Figure 4 with dots (note 
that these sets are not independent since each is a subset of 
the latter sets). A greater than 1 a deviation of a statistic 
rom the value expected for isotropy is indicated by the 
corresponding dot lying outside of the outer solid lines 
which represent the ± 1 a envelopes of the distributions of 
the statistics, including the effects of BATSE’s sky exposure. 
This occurs several times, but only by a small margin. Thus 
the depicted statistics are all consistent with isotropy. 

In order to test BATSE’s sky exposure map, Table 5 
presents the Earth-based statistics, <sin d) and <sin 2 S 
- t>, for the revised IB catalog, the 2B catalog, and the 


Number 


Data Set 
Name 


IB revised 
2B 

first 1005 


Size 

262 

585 

1005 


1 AnLh 5 

Results: Earth Coordinate System Tests 


Value 

0.015 

0.003 

0.016 


<sin <5> 

Deviation 3 Raw Deviation* 5 


-0.3 
— 0.9 
-0.5 


+ 0.4 
+ 0.1 
+ 0.9 


Value 

0.000 

0.024 

0.019 


<sin 2 <5 - 

Deviation 3 Raw Deviation 15 


-1.3 

- 0.1 

-0.7 


b Deviation, in <x, of the observed value from the value e'xpecS 


0.0 
+ 2.0 
+ 2.0 
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first 1005 GRBs. The largest deviation from the predictions 
based upon the sky exposure map is 1.3 a; the largest devi- 
ation for the first 1005 GRBs is only 0.7 <x. In contrast the 
value of <sin 2 d - for the 1005 GRBs is 2.0 a away from 
the mean that would be expected for a uniform sky expo- 
sure. This statistic is the most sensitive to BATSE’s nonum- 
form sky exposure and appears to have detected that effect 
at modest significance. These results are good indications 
that we are accurately predicting the effects of BATSEs 

nonuniform sky exposure. , n(V .. rD11[ 

Table 6 lists the subsets into which the set of 1005 GRBs 
have been split in order to test for possibly anisotropic 
subsets. Many of these sets have common members and are 
thus not independent. For each of these sets the Galactic 
and coordinate-system independent statistics have been cal- 
culated and are listed in Table 7. . . 

Because Table 7 lists the values of four statistics for 48 
subsets of the 1005 GRBs, it is likely that at least one value 
of a statistic will have a large deviation from the predictions 


of the null hypothesis of isotropy. The four statistics are not 
completely independent, just as the various subsets are not 
all independent. Ultimately, the degree of belief correspond- 
ing to a confidence level is a scientific judgement. In our 
judgement, for the entire data set an improbability of <0.01 
would be interesting and <0.001 would be convincing. 
Because of the many subsets examined, we would decrease 
these thresholds by at least a factor of 10 when considering 
a candidate anisotropy present only in a subset. 

Of the 48 data sets analyzed in Table 7, only one deviates 
from isotropy at more then the 99% confidence level, 
namely data set 46, the Medium GRBs of the revised IB 
catalog. These 51 GRBs are “Medium” according to the 
selection criteria of Quashnock & Lamb (1993a): log 10 
V < -0.8 and 465 <B< 1169, where B is a measure o! 
GRB brightness on the 1024 ms timescale and V compares 
the peak counts rates on the 64 and 1024 ms timescales 
(Lamb, Graziani, & Smith 1993). This data set deviates 
from isotropy by + 2.4 a for <cos 0) and by —2.6 a for 


Number Name 


1 

first 1005 

1005 

2 

IB revised 

262 

3 

2B 

585 

4 

2B less IB 

323 

5 

first 1005 less 2B 

420 

6 

2B less MAXBC 

485 

7 

first 1005 less MAXBC 

884 

8 

1st fifth 

201 

9 

2d fifth 

201 

10 

3d fifth 

201 

11 

4th fifth 

201 

12 

5th fifth 

201 

13 

r so and t g0 available 

551 

14 

r 50 short 

140 

15 

*50 l ° n S 

411 

16 

r qo short 

143 

17 

t 90 lon g 

408 

18 

flux/fluence available 

793 

19 

fluence < break 

718 

20 

fluence > break 

75 

21 

faintest fluence 

199 

22 

faint fluence 

198 

23 

bright fluence 

198 

24 

brightest fluence 

198 

25 

very bright fluence 

99 

26 

64 ms flux < break 

734 

27 

64 ms flux > break 

59 

28 

faintest 64 ms flux 

199 

29 

faint 64 ms flux 

198 

30 

bright 64 ms flux 

198 

31 

brightest 64 ms flux 

198 

32 

very bright 64 ms flux 

99 

33 

1024 ms flux < break 

770 

34 .... 

1024 ms flux > break 

23 

35 .... 

faintest 1024 ms flux 

199 

36 .... 

faint 1024 ms flux 

198 

37 .... 

bright 1024 ms flux 

198 

38 .... 

brightest 1024 ms flux 

198 

39 .... 

very bright 1024 ms flux 

99 

40 .... 

soft events 

261 

41 .... 

intermediate hardness 

261 

42 .... 

hard events 

261 

43 .... 

C /C min available 

601 

44 .... 

64 ms triggers 

418 

45 .... 

1024 ms triggers 

504 

46 ... 

IB medium 

51 

47 .... 

first 1005 less IB medium 

76 

48 ... 

first 1005 medium 

127 


TABLE 6 
The Data Sets 


Description 


triggers 105 3005, 1991 April 19-1994 May 31. Post-2B locations are preliminary 

the revised IB catalog, triggers 105-1466 (Meegan et al. 1994a) 

the 2B catalog, triggers 105-2230. (Meegan et al. 1994a) 

triggers 1467-2230 

triggers 2232-3005 

triggers 105-2230 with MAXBC-based locations removed 

triggers 105-3005 with MAXBC-based locations removed 

triggers 105-1204 

triggers 1211-1723 

triggers 1726-2280 

triggers 2281-2606 

triggers 2608-3005 

t 50 and t 9Q have been calculated 

ho < 0 8 s 
f 50 > 0.8 s 
t 90 < 2.0 s 

flux and fluence numbers available; a subset of first 1005 ^ 

inhomogeneous: 50-300 keV fluence < break at 1.0E-5 erg cm^ 
homogeneous: 50-300 keV fluence > break at 1.0E-5 erg cm 
faintest quarter by 50-300 keV fluence: 1.1E-8 to 2.7E-7 erg cm 
2nd faintest quarter by 50-300 keV fluence : 2.7E - 7 to 9.5 E - 7 erg cm ^ 

2nd brightest quarter by 50-300 keV fluence: 9.6E-7 to 3.2E-6 erg cm 
brightest quarter by 50-300 keV fluence: 3.3E — 6 to 1.4E— -4 erg ^ 
brightest eighth by 50-300 keV fluence : 7.8E - 6 to 1 ,4E - 4 erg cm _ 
inhomogeneous: 50-300 keV 64 ms peak flux < break at 10. y s cm 
homogeneous : 50-300 keV 64 ms peak flux > break at 10. y s cm 
faintest quarter by 50-300 keV 64 ms peak flux: 0.32 to 0.87 y s ctn ^ 
2nd faintest quarter by 50-300 keV 64 ms peak flux : 0.88 to 1 .45 y s an 
2nd brightest quarter by 50-300 keV 64 ms peak flux: 1-45 to 3.00 y s ^ 
brightest quarter by 50-300 keV 64 ms peak flux: 3.00 to 168. y s cm ^ 
brightest eighth by 50-300 keV 64 ms peak flux: 6.18 to 168 y s cm ^ 
inhomogeneous : 50-300 keV 1024 ms peak flux < break at 1 3. ys^ cm 
homogeneous: 50-300 keV 1024 ms peak flux > break at 13. y s cm ^ 
faintest quarter by 50-300 keV 1024 ms peak flux: 0.05 to 0.40 y s cm 2 
2nd faintest quarter by 50-300 keV 1024 ms peak flux : 0.40 to 0.69 y s cm 
2nd brightest quarter by 50-300 keV 1024 ms peak flux. 0.69 to 1.60 y s cm 
brightest quarter by 50-300 keV 1024 ms peak flux: 1.60 to 39.2 y * C » J 
brightest eighth by 50-300 keV 1024 ms peak flux: 3.45 to 39.2 y s cm 
50-100 keV/ 100-300 keV fluence < 0.2611 
50-100 keV/ 100-300 keV fluence between 0.2612 & 0.4411 
50-100 keV/ 100-300 keV fluence > 0.4418 
data available to calculate C m JC min 

events that exceeded trigger threshold on 64 ms timescale, subset of 43 
events that exceeded trigger threshold on 1024 ms timescale subset of 43 
triggers < 1466 satisfying Medium definition of Quashnock & Lamb (1993a) 
triggers 1467-3005 satisfying Medium definition of Quashnock & (1993a) 

triggers <3005 satisfying Medium definition of Quashnock & Lamb (1993a) 
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TABLE 7 


Results: Galactic and Coordinate System Independent Tests 



Data Set 


<cos 9} 

<sin 2 

b-i) 


it 


£ 

Number 

Name 

Size 

Value* 

Deviation 15 

Value* 

Deviation* 5 

Value* 

Deviation* 5 

Value* 

Deviation* 5 

1 

first 1005 

1005 

0.004 

+ 0.9 

-0.008 

1 

© 

1.6 

-0.9 

6.2 

-1.0 

2 

IB revised 

262 

0.026 

+ 1.1 

-0.025 

-1.1 

1.5 

-0.7 

6.0 

-0.2 

3 

2B 

585 

0.022 

+ 1.5 

-0.013 

-0.7 

1.5 

-0.8 

6.5 

-0.6 

4 

2B less IB 

323 

0.019 

+ 1.0 

-0.004 

+ 0.1 

0.4 

-1.1 

7.0 

-0.1 

5 

first 1005 less 2B 

420 

-0.022 

-0.3 

-0.002 

+ 0.2 

1.6 

-0.7 

5.2 

-0.6 

6 

2B less MAXBC 

485 

0.027 

+ 1.5 

-0.025 

-1.5 

2.1 

-0.6 

5.6 

-0.6 

7 

first 1005 less MAXBC 

884 

0.005 

+ 1.0 

-0.013 

-0.8 

2.6 

-0.6 

5.6 

-1.0 

8 

1st fifth 

201 

0.037 

+ 1.2 

-0.033 

-1.4 

1.2 

-0.8 

7.4 

+ 0.2 

9 

2d fifth 

201 

0.052 

+ 1.6 

0.023 

+ 1.3 

4.8 

+ 0.5 

14.1 

+ 1.9 

10 

3d fifth 

201 

-0.024 

-0.3 

-0.031 

-1.2 

1.5 

-0.7 

6.5 

0.0 

11 

4th fifth 

201 

0.001 

+ 0.3 

-0.014 

-0.4 

0.1 

-1.2 

4.5 

-0.5 

12 

5th fifth 

201 

-0.048 

-0.8 

0.013 

+ 0.8 

4.8 

+ 0.5 

2.9 

-0.9 

13 

t 50 and t 90 available 

551 

0.017 

+ 1.2 

-0.021 

-1.3 

3.7 

-0.1 

4.9 

-0.8 

14 

f 50 short 

140 

-0.010 

+ 0.1 

-0.025 

-0.8 

0.5 

-1.0 

3.2 

-0.8 

15 

Go l° n g 

411 

0.026 

+ 1.4 

-0.020 

-1.0 

4.2 

+ 0.1 

2.7 

-1.1 

16 

t 90 short 

143 

-0.007 

+ 0.1 

-0.018 

-0.5 

0.3 

-1.1 

2.8 

-0.9 

17 

Go lon g 

408 

0.026 

+ 1.4 

-0.022 

-1.2 

4.6 

+ 0.2 

3.7 

-0.9 

18 

flux/fluence available 

793 

0.011 

+ 1.2 

-0.012 

-0.7 

1.0 

-1.0 

2.9 

-1.4 

19 

fluence < break 

718 

0.015 

+ 1.3 

-0.015 

-0.9 

2.0 

-0.7 

3.0 

-1.3 

20 

fluence > break 

75 

-0.026 

-0.2 

0.012 

+ 0.5 

2.0 

-0.4 

4.7 

-0.2 

21 

faintest fluence 

199 

0.015 

+ 0.7 

0.003 

+ 0.4 

2.6 

-0.3 

8.3 

+ 0.4 

22 

faint fluence 

198 

-0.016 

-0.1 

-0.015 

-0.5 

1.6 

-0.6 

5.0 

-0.4 

23 

bright fluence 

198 

0.018 

+ 0.8 

-0.034 

-1.4 

1.3 

-0.7 

6.2 

-0.1 

24 

brightest fluence 

198 

0.026 

+ 1.0 

-0.004 

+ 0.1 

3.3 

0.0 

0.8 

-1.4 

25 

very brightest fluence 

99 

0.040 

+ 0.9 

0.006 

+ 0.4 

3.0 

-0.1 

4.0 

-0.4 

26 

64 ms flux < break 

734 

0.000 

+ 0.6 

-0.012 

-0.6 

1.1 

-1.0 

2.6 

-1.4 

27 

64 ms flux > break 

59 

0.146 

+ 2.1 

-0.018 

-0.4 

6.6 

+ 1.4 

1.3 

-1.2 

28 

faintest 64 ms flux 

199 

0.023 

+ 0.9 

0.001 

+ 0.3 

2.2 

-0.4 

2.8 

-0.9 

29 

faint 64 ms flux 

198 

0.048 

+ 1.5 

-0.012 

-0.3 

1.4 

-0.7 

5.8 

-0.2 

30 

bright 64 ms flux 

198 

-0.052 

-1.0 

-0.007 

-0.1 

7.1 

+ 1.3 

5.4 

-0.3 

31 

brightest 64 ms flux 

198 

0.026 

+ 1.0 

-0.031 

-1.2 

4.3 

+ 0.3 

3.0 

-0.9 

32 

very bright 64 ms flux 

99 

0.068 

+ 1.4 

-0.054 

-1.7 

4.4 

+ 0.4 

5.1 

-0.2 

33 

1024 ms flux < break 

770 

0.006 

+ 0.9 

-0.014 

-0.8 

1.2 

-1.0 

3.0 

- 1.3 

34 

1024 ms flux > break 

23 

0.166 

+ 1.5 

0.031 

+ 0.6 

2.8 

-0.1 

1.2 

-12 

35 

faintest 1024 ms flux 

199 

-0.039 

-0.6 

0.014 

+ 0.9 

1.0 

-0.8 

7.5 

+ 0.2 

36 

faint 1024 ms flux 

198 

0.061 

+ 1.8 

-0.021 

-0.8 

5.4 

+ 0.7 

7.3 

+ 0.2 

37 

bright 1024 ms flux 

198 

0.038 

+ 1.2 

-0.002 

+ 0.1 

1.4 

-0.7 

6.5 

0.0 

38 

brightest 1024 ms flux 

198 

-0.016 

-0.1 

-0.041 

-1.7 

2.5 

-0.3 

5.1 

-0.4 

39 

very bright 1024 ms flux 

99 

0.068 

+ 1.4 

-0.030 

-0.8 

1.4 

-0.6 

4.9 

-0.2 

40 

soft events 

261 

0.008 

+ 0.6 

-0.024 

-1.0 

2.0 

-0.5 

5.9 

-0.3 

41 

intermediate hardness 

261 

0.046 

+ 1.7 

-0.019 

-0.8 

2.0 

-0.5 

2.5 

-1.0 

42 

hard events 

261 

-0.026 

-0.4 

0.002 

+ 0.4 

1.7 

-0.6 

0.3 

-1.6 

43 

C m „/c mi „ available 

601 

-0.021 

-0.4 

-0.026 

-1.7 

1.0 

-1.0 

8.9 

- 0.1 

44 

64 ms triggers 

418 

-0.017 

-0.1 

-0.033 

-1.9 

2.4 

-0.5 

8.4 

0.0 

45 

1024 ms triggers 

504 

-0.017 

-0.1 

-0.026 

-1.6 

0.5 

-1.1 

6.4 

-0.5 

46 

IB medium 

51 

0.184 

+ 2.4 

-0.111 

-2.6 

5.3 

+ 0.9 

7.3 

+ 0.6 

47 

first 1005 less IB medium 

76 

-0.127 

-1.7 

-0.014 

-0.3 

4.6 

+ 0.6 

8.2 

+ 0.8 

48 

first 1005 medium 

127 

-0.002 

+ 0.2 

-0.053 

-1.8 

0.2 

-1.1 

9.1 

+ 0.8 


* Observed value uncorrected for BATSE’s sky exposure. 

b Deviation, in cr, of the observed value from value calculated from BATSE's sky exposure map assuming isotropy. 


<sin 2 b — ^>. These deviations are in the direction of a con- 
centration toward the Galactic center and the Galactic 
plane. The data set of Quashnock & Lamb (1993a) was 
selected from the IB catalog and had four additional GRBs; 
the C max /C min -values of these four GRBs on untriggered 
timescales have been withdrawn by the BATSE team as 
being erroneous and these four GRBs are, therefore, not 
included in the Medium bursts of the revised IB catalog. 
The 55 GRB data set of Quashnock & Lamb (1993a) devi- 
ates from the predictions of isotropy and BATSE’s sky 
exposure map by +3.2 o for <cos 6 ) and by —2.9 a for 
<sin 2 b - !>. The 51 Medium GRBs of the revised IB 
catalog are shown with filled circles on Figure 7 while the 
additional four from the original 1 B catalog are shown with 
open circles. 

The dipole and quadrupole moments of the Quashnock 


and Lamb 55 GRB data set each have significances of 
roughly 10“ 3 . Quashnock and Lamb point out that these 
two statistics are partially independent and find from 
Monte Carlo simulations that the probability of seeing both 
deviations or larger is 1.9 x 10 6 . Because they optimized 
the brightness B selection criterion to obtain this result, 
they estimate the probability as 1.1 x 10“ 4 . They inter- 
preted their data set (filled and open circles in Fig. 7) as 
evidence for spiral arms and concluded that all gamma-ray 
bursts are Galactic in origin (Quashnock & Lamb 1993a). 

The estimation of the correction factor to apply to signifi- 
cances obtained from retrospective analysis of data is diffi- 
cult. How many ways were the data subdivided? How many 
statistical approaches were tried? In this case, we should 
have some additional factor due to the V selection criterion 
and to the idea of using a joint probability from <cos 9 > 
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Fig. 7. — GRBs in the Medium category of Quashnock & Lamb 
(1993a). The locations marked with circles are their data set. The open 
circles mark the four GRBs of the IB catalog, which can no longer be 
classified according to the Medium criterion because the C m „/C min -values 
for these GRBs on nontriggered timescales have been withdrawn by the 
BATSE team as erroneous. The filled circles are the Medium GRBs of the 
revised IB catalog (data set 46) and the stars are the post- IB Medium 
GRBs (data set 47). The map is a Hammer-AitolT projection in Galactic 
coordinates. 


and <sin 2 b — ^>. While retrospective analysis has poten- 
tial pitfalls, it is necessary because of our limited under- 
standing of gamma-ray bursts. 

The best solution is to collect more data in order to make 
an independent test of a hypothesis formed by retrospective 
analysis. Using exactly the same selection criteria as Quash- 
nock and Lamb (Graziani 1993), we find 76 additional 
Medium GRBs in the first 1005 GRBs — these form our data 
set 47 and are shown with stars on Figure 7. The rate of 
discovering Medium GRBs has decreased due to data gaps 
caused by the failure of the CGRO tape recorders, which 
prevent our determining the B- and V - values for many 
GRBs. Even a small data gap, which does not hinder GRB 
localization, prevents determination of the peak count rate 
because the peak rate may have occurred during the data 
gap. The data gaps should not have any correlation with 
intrinsic properties of gamma-ray bursts and thus should 
not bias the selection. The new Medium GRBs are consis- 
tent with isotropy and in fact the —1.7 a deviation of 
<cos 0) indicates a deficiency of GRBs toward the Galactic 
center. The evidence for a Galactic origin of GRBs is not 
confirmed by a larger, independent sample. The entire data 
set of 51 + 76 = 127 Medium GRBs, data set 48, is also 
consistent with isotropy: <cos 0 ) deviates by +0.2 a and 
<sin 2 b — j) by — 1.8 a. 

What can we conclude from the contradictory results of 
data sets 46 and 47? If one takes data set 46 to be evidence 
of a true anisotropy and thus the values of <cos 0 > and 
<sin 2 b — 3> observed for this set to be measurements of 
true anisotropic moments, then one has the difficulty of 
explaining the 3.0 cr and 1.8 a deviations of the moments 
observed for data set 47 from these values. Conversely, if 
one assumes the GRBs to be isotropic, then one has the 
difficulty of explaining the + 2.4 a and the — 2.6 o of devi- 
ations of <cos 0) and <sin 2 b — j> for data set 46 from the 
values expected for isotropy. Because data set 46 was identi- 
fied by retrospective analysis and its result is not confirmed 
by a larger, independent sample, we judge that the evidence 
against the null hypothesis of isotropy is weak. 

We use the statistics of Horack et al. (1994) to search for a 
heliocentric pattern in the first 1005 GRBs. The dipole 
moment toward the Sun, <cos (/>), is observed to be 0.006, 
which deviates by +0.3 a from the value of 0.000 predicted 


by isotropy and BATSE’s sky exposure. Because the direc- 
tion to the Sun changes with time, the sky exposure bias is 
calculated by averaging the time-dependent bias specified in 
Table 3 over 3.1 yr. The quadrupole moment in the ecliptic 
plane, <sin 2 /? — 3), is observed to be 0.022, a +0.2 a devi- 
ation from the value of 0.020 predicted by isotropy and 
BATSE’s sky exposure. No evidence for a heliocentric 
origin is found. 

We have searched for an excess or deficiency of GRBs 
near the Magellanic Clouds, the Virgo cluster, M31, the 
Galactic center, and a Centauri. In all 48 data sets there is 
no instance in which the number of GRBs within 20° of 
these objects is less probable than 1%. 

In summary, we have examined 48 subsets of the first 
1005 gamma-ray bursts observed by BATSE using four 
dipole and quadrupole tests. The entire set of 1005 GRBs 
has been tested with four additional dipole and quadrupole 
statistics. Only one set, 46, appears anisotropic, however, a 
larger, independent set selected with the same criteria, 47, 
fails to confirm the anisotropy so we judge the data to be 
consistent with isotropy. The largest anisotropy in any set 
other than 46 is the 2.1 a deviation of <cos 0) for data set 
27. Because of the large number of partially independent 
sets and statistics examined, a deviation of this magnitude is 
expected. The gamma-ray bursts locations are consistent 
with large-scale isotropy. 

5. MODEL COMPARISONS AND LIMITS 

While the statistics show that gamma-ray burst locations 
are consistent with isotropy, they do not prove that gamma- 
ray bursts are isotropic— one could imagine that a small 
anisotropy exists that cannot be detected in a sample of 
1005 locations. This raises the question of how large an 
anisotropy could exist and not have been detected. We 
address this question in two ways : First, we discuss below 
and summarize in Table 8 the comparisons of the moments 
observed by BATSE with the moments of published Galac- 
tic models. Secondly, we use Monte Carlo simulations to 
determine, for a variety of models, which anisotropies could 
have been detected. These upper limits on anisotropies are 
discussed below and summarized in Figure 8. 

For both the model comparisons (Table 8) and model 
parameter limits (Fig. 8) we have assumed that the number 
of sources is the same as the number of observed bursts. If 
some of the sources repeat so that there are fewer sources 
than observed bursts, the finite sample fluctuations will be 
larger than we have calculated, which will reduce our ability 
to detect anisotropies (Quashnock 1995). We have placed a 
20% upper limit on the repeater fraction at the 99% con- 
fidence level (Meegan et al. 1995), which implies that there 
are at least nine sources for each 10 observed bursts. In this 
case the finite sample fluctuations would be 10% larger than 
we have calculated. 

5.1. Model Comparison Methods 

Below and in Table 8 we compare the moments of 
published Galactic models with the observations. Because 
the goal of these comparisons is to confront the models with 
BATSE’s observations, the comparisons are made between 
the observed moments corrected for the sky exposure bias 
and the model moments. These published models were 
developed based upon earlier and smaller data sets — some 
could have improved agreement with the current observ- 
ations by reoptimizing their parameters. In other cases, full 
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TABLE 8 


Post-BATSE Galactic Models 


Model 


Statistic 

Model Value 

Deviation 1 

Eichler & Silk (1992) 
Hartmann (1992) . 


<cos 0 > 

/ 0 ;„2 l i \ 

0.05 

1.8 

Li & Dermer (1992) . 


\Sin o — -^2 

— 0.05 

5.0 

Lingenfelter & Higdon (1992) 

Fabian & Podsiadlowski (1993) 

Smith & Lamb (1993): disk/Gaussian halo 

Smith & Lamb (1993): dark matter halo/disk 

Higdon & Lingenfelter (1994), R core = 7.5 kpc, 25% disk... 

Higdon & Lingenfelter (1994), R eott = 15 kpc, 20% disk 

Higdon & Lingenfelter (1994), R cort = 30 kpc, 8% disk 

Li et al. (1994) c 

\COS V) 

<cos 0} 
<sin 2 b - 
<cos 0 LMC >" 
<sin 2 b — 
<cos 0 > 
<cos 6} 
<cos 0> 
<cos 0 > 

0.048 

0.08 

-0.06 

0.038 

-0.027 

0.057 

0.088 

0.073 

0.060 

1.7 

3.5 
6.1 

1.9 

2.6 
2.2 

3.9 
3.1 
2.4 

Podsiadlowski et al. (1995), Fig. 5 a 


\Sin D — 
<cos 2 6 — ^> d 

/ PA c \ 

- 0.084 
0.073 
0.043 
-0.019 
0.054 
-0.024 
0.056 

2.0 

3.1 

Podsiadlowski et al. (1995), Fig. 5b . 
Smith (1994b) .... 


U / 

<sin 2 b — 
<cos 6 ) 
<sin 2 b - ■£> 
<cos 0> 

1.4 

1.7 

2.0 

2.2 

2.1 


I ncl u^es rarrectio n for £ ^ ° bSerVed ^ 1005 bUrS * S " burs,s for Li e * al >**)• 

sky'exjwsure 6^0^“ '° ^ ^ a ° ud = * he observed value is -0.020, and the predicted 

set S w7h99 burT'" ^ ** Wi ' h 1024 mS ^ fluX greater than 145 1 cm ‘ 2 , and the comparison data 

d The observed value of this quadruple moment is -0.023, and the sky exposure predicted value is -0.004. 


reoptimization might increase the deviations in the 
moments because of the tighter constraints of the brightness 
distribution data. Some of the models are merely examples 
of good fits to older data, rather than best-fits to that data. 
These models might be improved to have reduced moments, 
possibly at the expense of less “reasonable” parameters’ 
The uncertainties in the theoretical models are difficult to 
quantify. In any case, reoptimizing the models is beyond the 
scope of this paper, so we compare the observations with 
the published moments without any corrections. We 
convert the differences between the observed and published 
model moments to deviations in units of a, using for a the 
finite sample fluctuations of the observations, ignoring the 
possibility of uncertainties in the models. 

The values of a are calculated as the sample fluctuations 
expected assuming isotropy (Table 3)— a more accurate 
procedure would be to calculate the values of a from the 
anisotropic models (Li 1995). For models with moments so 
small that ~ 1000 bursts are needed to detect the moment, 
the difference between the model and isotropy c-values are 
typically <1%. For the model of Li et al. (1994), which 
predicts a larger moment in a subset of the data, the differ- 
ence somewhat exceeds 10% (Li 1995). Since the differences 
between the isotropy and model er-values are small and the 
model ff-values generally have not been published, we cal- 
culate the <T-values assuming isotropy. 

In § 3.2 we found that smearing the locations by the 
angles £ listed in Table 2 reduced the dipole moment 
<cos^0> by the factor 0.86 and the quadrupole moment 
<sin b — j) by 0.68. These factors considerably exaggerate 
the effects of smearing since the angles I are chosen to be at 
least twice ff (ot for the least well-located GRB in each of four 
categories (Table 2). Simulations using more accurate 
smearing angles show the smearing factors to be only a few 
percent below 1, and we have, therefore, neglected this effect 
in our model comparisons (Table 8). 


5.2. Model-Limit Determination Methods 

Monte Carlo simulations are used to determine for 
various models, as a function of the number of bursts N B , 
the parameter values that should create detectable aniso- 
tropies. The goal of the simulations is to determine what 
anisotropies should be detectable and, therefore, the upper 
limits are not based upon the moments observed bv 
BATSE. y 

In order to calculate the ability of a statistic to detect an 
anisotropy, the “power” of a statistic, it is necessary to 
assume a particular anisotropic model (Eadie et al. 1971 ; 
Martin 1971). For example, the statistic <cos 0> is highly 
sensitive to a dipole moment toward the Galactic center but 
completely insensitive to a dipole moment toward the north 
Galactic pole. For each model we determine the detect- 
ability of anisotropies for three confidence levels, 99% 
99.9%, and 99.99%, which are equivalent to 2.3 <r, 3.1 <x, and 
3.7 ff for Gaussian-distributed statistics. The critical values 
of the statistics are determined from the distributions of the 
statistics under isotropy so that the chance probabilities of 
mistakenly deciding that an isotropic distribution is aniso- 
tropic is 1%, 0.1%, and 0.01% (see the vertical lines on Fig. 
5). Choosing a high confidence level reduces the chance of 
falsely discovering an anisotropy, a mistake known as a 
type I error, at the expense of reducing the ability to dis- 
cover a true anisotropy (Eadie et al. 1971 ; Martin 1971). We 
show the three confidence levels to allow readers to choose 
based upon their judgement and considerations such as 
whether the entire data set or a subset is being tested. 

Just as an isotropic distribution may by chance look 
rather anisotropic, so may an anisotropic distribution 
appear rather isotropic. This is seen in Figure 5 in the areas 
of the anisotropic histograms that extend near the center of 
the isotropic histograms. A Type II error is failing to dis- 
cover an anisotropy when it actually exists and the power of 



Fig. 8a 



Fig. 8 b 

F,C. S .-L- S o„,h,»^ot r «,,^^ 

discovered, i.e., is found in 50 /« of the simulations. I ne a assumed visible ( 6 ) Watson distribution. All GRBs are assumed visible, (c) 

P f ° Und '".d 5 :! The^GR^are'distribut^expOTerU^Hy^ndistam^fromThe^l^ne^The distribution is parameterized by the observing depth, which is the 
Exponential disk. The GRBs are di istnbuted expone y The GRB denslt tg umform m an in fi n itely long 

model. A mixture of /, ni „ locations drawn from a Watson distribution with k — with the remainder isotropic. 
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Fig. 8< 



Fig. 8 d 


a statistic is one minus the probability of making a Type II 
error, i.e., the probability of detecting the anisotropy (Eadie 
et al. 1971; Martin 1971). We choose to determine the 
detectability of anisotropies for two powers: for BATSE 
typically and “nearly always” discovering an anisot- 
ropy, by which we mean finding the anisotropy at or past 
the stated confidence level in 50% and 95%, respectively, of 
the simulations. 

Figure 8 shows as a function of the number of gamma-ray 


bursts N b the values of the anisotropy parameters such that 
either 50 A or 95% of the simulations find the anisotropy at 
one of our confidence levels. The detectabilities are deter- 
mined as a function of N B for N„ = 100 to 2000. This allows 
one to determine what anisotropy could be found in a 
subset of the current 1005 GRBs or will be detectable in the 
future. Each panel of Figure 8, thus, has six curves: three 
confidence levels at two different powers. The “typical” 
detection curves are solid, while the “ nearly always ” curves 
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Fig. 8e 



Fig. 8/ 


are dashed. The y-axes are arranged so that more aniso- 
tropic models are at the top of the graph, hence, the top 
curves are for the 99.99% confidence level and the bottom 

curves are for the 99% confidence level. 

We note that in all cases detecting an anisotropy nearly 
always” at the 99% confidence level ( lowest dashed curves) 
requires a larger anisotropy than detecting the anisotropy 
“typically” at the 99.99% confidence level ( highest solid 
curve). This is an indication of how conservative in the sense 


of minimizing BATSE’s ability to find anisotropies the 

“ nearly always ” requirement is. 

The detectability of each model is determined for the 
most sensitive Galactic statistic for that model ( left panels) 
and for the most sensitive coordinate-system independent 
statistic ( right panels). This allows one to determine what 
Galactic anisotropies and what anisotropies in less expected 
directions, e.g., M31, Virgo cluster, Magellanic Clouds, 
BATSE could detect. With one exception discussed below 
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(§ 5.4), the Galactic-based statistic are more sensitive for 

G^ZZrLZT ^ n,odels - which are cons,ruc " d 

The location errors slightly “isotropize” the observ- 
ations. The Monte Carlo simulations used to determine the 
detectability of the model anisotropies (Fig. 8) include a 
smearing of the simulated locations by an angle a = 14 °. 
The value 14 is the root mean square sum of the 13° sta- 
tishca 1 and 4 systematic error typical for BATSE’s faint 
CRBs and is, therefore, conservative in the sense of under- 
estimating our ability to find anisotropies. 

5.3. Fisher and Watson Models 
We use the Fisher and Watson distributions as “generic ” 
dipo e and quadruple distributions to derive some illustra- 
tive limits. We assume that all locations generated in the 
simulations of these distributions are “observed” Figures 
8a and 8 b show the detectable values of the concentration 
parameter * : for these two models. Values of k may be con- 
verted into dipole and quadruple moments and into angles 
containing one-half of the probability using Figure 2 
The lowest curve of the left box of Figure 8a shows that 
or 1000 GRBs BATSE would “typically,” i.e., 50% of the 
tirae,dc(ect a Fisher distribution with k > 0.13 at the 99% 
confidence level or higher; a simulated example of this dis- 

Jrr,TftT r , !k Sh ^ Wri -u n FigUre la ' The di P° le moment 
<cos 0> of this distribution is 0.043, implying an expected 

y ificance for lOOO GRBs of 2.4 a. The Fisher distribution 
^ has one-half of the probability within 0.,, = 

86 of the Galactic center compared to 90° for isotropy 
A very conservative upper limit obtained from the middle 
dashed curve of the left box of Figure 8a is that BATSE 
would with 1000 locations “nearly always,” i.e., 95% of the 
de ‘ e ^4 Fish u er distribution about the Galactic center 
th k > 0.27 at the 99.9% confidence level or better. This 
nio tnb , Ut '£ n one-half of the probability within 0 U2 = 

82 of the Galactic center and an expected moment 


r»l 8> .Z° 090 ' implying a typical detection for 1000 
OKBs at the very high level of 4.9 a. 

* S ™ ll £ result s/or finding a quadrupole moment created 
by the Watson distribution are given in Figure 8 b. A value 
ot k - -0.28 would typically be discovered in 1000 loca- 

FM ,‘m tk”? 7-P ,eS V at ,h ' W% level 

(see Fig. 16). This distribution has one-half of the probabil- 
ity within J>i /2 - 28° of the Galactic plane versus 30° for 
isotropy The fact that a plane concentration this weak is 
detecdabte but not detected shows how isotropic the GRBs 
are. The middle dashed curve of the left box of Figure 86 
ratcc 6 ^ c °" servative u PP er limits, e.g., for 1000 GRBs 

Ipv J SE «/ W ? U d !? Car ^ a ' Ways find at the 99 - 9% confidence 
level a Watson distribution with k = - 0.57. This distribu- 

tion has one-half of the probability within 26° of the Galac 
tic plane. 

The right box of Figure 8 a gives results for detecting a 
concentrahon about any point on the sky: a value of 
D , •? ^? ud ** ‘yprcally detected in 1000 GRBs by the 
Ray eigh- Watson dipole statistic tT at the 99% confidence 
evel. The corresponding concentration is having one-half of 
the probability within 85° of any point. Similarly, the right 
box of Figure 86 gives results for detecting a concentratfon 
about any plane on the sky: a value of k < -0.4 would be 
yp.cally detected in 1000 GRBs by the Bingham quadrupo- 
le statistic 36 at the 99% confidence level. The correspoL- 

:^rr° n ls , havin 8 one -half of the probability 
within 27 of some plane. These results constrain sugges- 
uons that gamma-ray bursts originate from the Magellanic 
Clouds (Maoz 1993a; Fabian & Podsiadlowski 1993) or 
rom the Milky Way and M31 (Gurevich et al. 1994) or 
from any other unexpected direction on the sky. 

5.4. Disk and Spiral Arm Models 
The Watson distribution discussed above serves as a non- 
physical model of a disk distribution. We also consider a 
somewhat more physical model: the GRBs are assumed to 
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be standard candles distributed exponentially in distance 
from the Galactic plane and the solar system is assumed o 
lie in the plane. Note that the vertical profile of the thin 
stellar disk of the Milky Way is adequately modeledbyan 
exponential (Gilmore, King, & van der Kruit 1 )■ 

appearance of this distribution is parameterized by the ratio 
observing depth, which is the maximum distance to whl ™ 
GRBs can be seen, divided by the exponential scale heig _ 

For 1000 GRBs, even using the very conservative limit o 
nearly always detecting the anisotropy at the 99. /o con 
fidence level, the middle dashed curve of Figure 8c ef 
shows that BATSE cannot be viewing past 0.8 exponents 
scale heights. In this model BATSE is viewing a spherica 
region no larger than 0.8 scale heights and the ratio of the 
lowest density of GRBs viewed over the highest density is 
exp ( — 0.8) = 0.45, which is reached only at the north and 
south Galactic poles of the spherical region. Since the faint 
GRBs are actually deficient by a factor of about 6, t is 
model clearly does not agree with the inhomogeneous 
radial distribution observed by BATSE. 

As a very simple spiral arm model we consider an infi- 
nitely long cylinder containing a uniform density ot stan- 
dard candle sources, with the Sun located on the cylinder 
axis. This distribution is parameterized by the ratio observ- 
ing depth divided by arm radius. Using the very conserva- 
tive upper limit of nearly always detecting the anisotropy at 
the 99 9% confidence level, the middle dashed curve o 
Figure 8<f right shows that the observing depth cannot 
exceed 1.2 arm radii. The fractional volume of the sphere in 
which GRBs could be seen occupied by the uniform density 
cylinder is 83%. Again, this clearly does not match the 
major deficiency of faint GRBs observed by BATS 

Only in the spiral arm model is the coordinate-system 
independent statistic more sensitive than the Galactic-based 
statistic. This is because <sin 2 b - j> is only detecting the 
concentration the arm makes in the plane, while the 
Bingham statistic cJS is also detecting the concentration ol 
GRBs at the two poles of the cylinder axis. We could con- 
struct a special Galactic-based statistic to look for this latter 
effect, but the limits obtained with iM amply suffice. 

In both the disk and spiral arm models one would expect 
the faintest sources to show the largest anisotropy, but no 
anisotropy is detected in these GRBs, e.g., data sets 21, 28, 

35 Jt should be clear that a disk or spiral arm model requires 
a strong anisotropy because of the major deficiency of faint 
GRBs seen by BATSE. For example, in a disk model, he 
many faint GRBs that are “missing” above and below the 
plane in order to create the observed inhomogeneity would 
necessarily produce a concentration in the plane and a large 
quadrupole moment. This was shown to be the case for 
several disk models by Paczynski (1990,1991b) and has 
been confirmed for a variety of models by Mao & Paczynski 
(1992b) and Hakkila et al. (1994a). The inability of spira 
arm models to fit the BATSE data has been shown by 
Hakkila et al. (1994a) and by Smith (1994a). A model with 
more than one spiral arm visible would be subject to the 
constraints against an excess in the plane characteristic o 
disk models. 

5.5. Geometric Halo Models 

The basic requirements on a halo model are (1) the length 
scale of the distribution should be S> the solar galactocen- 
tric distance R 0 = 8.5 kpc, so that the distribution will 


appear sufficiently isotropic; (2) observing depth (the 
maximum distance seen by BATSE) should be S> the length 
scale of the distribution so that there will be a sufficient 
deficiency of faint GRBs, i.e., inhomogeneity; and (3) the 
observing depth should be <350 kpc, half the distance o 
M31, so that the model will not predict an excess of GRBs 
from M3 1. Since the focus of this paper is large-scale isot- 
ropy, we will primarily be concerned with the first require- 
ment- that our offset from the center of the galaxy be small 
compared to the size of the observed GRB distribution so 
that the dipole moment not be too large. Generally the 
constraint imposed by the quadrupole moment is less hmit- 

10 Any spherically symmetric halo model may be thought of 
as the sum of a series of shells. Hartmann et al. (1994) have 
shown that a shell of galactocentric radius R has a dipole 
moment 


2 Rq 

<cos 0} = J 


( 1 ) 


Thus, a shell of radius R = 1 10 kpc has a dipole moment 2o 
above the value observed for 1005 GRBs (corrected for the 
sky exposure bias). Any GRBs interior to this radius will 

have to be balanced by GRBs exterior to this radius. 

As our simplest halo model we consider a sphere centered 
on the Galactic center containing a uniform density ot 
sources, all of which are observable by BATSE. Wasserman 
(1992) has considered a similar model. The lowest so l 
curve of the left panel of Figure 8e shows that the radius 
r of this sphere must be greater than 194 kpc to avoid 
BATSE typically detecting the anisotropy in 1005 GRBs at 

the 99% confidence level. . , . . , , rm 

Our next model is based upon a simple functional form 
for the dark matter halo used to explain the rotation curves 
of galaxies including the Milky Way : density 


p(R) oc - 


+ (R/KcoJ 2 


( 2 ) 


The R -2 dependence yields the observed flat rotation 
curves, and the core radius prevents an excess contribution 
at the center and prevents the function from diverging. In 
contrast, the distribution of the luminous halo of our 
Galaxy is at least as steep as R (Gilmore et al. 1990, 
Djorgovski & Meylan 1994), which is more anisotropic 
than the R 2 distribution that we use. Hakkila et al (1994b) 
found using 452 GRBs that exponents from 1.5 to 4.0 were 
acceptable. Paczynski (1991b) and Mao & Paczyns i 
(1992b) used this distribution and concluded that the e 
would be an observable deficiency of faint sources and an 
undetectable dipole moment only if R core > 2 R 0 - 17 kpc_ 
Brainerd (1992) found that a somewhat smaller value o! 

R was permitted if GRBs are not standard candles 

Figure 8/ presents our results for a dark matter halo of 
standard candle sources observed to 300 kpc. This observ- 
ing depth is near the maximum permitted without seeing 
sources from M31 and is near the best value found by 
Hakkila et al. (1994a), considering both isotropy and mho- 
mogeneity. A smaller observing depth for a fixed number ot 
observed GRBs would lead to a larger dipole moment. We 
find that a core radius of R core = 23 kpc would be typically 
found in 1005 locations by BATSE at the 99% confidence 
level. Hakkila et al. (1994a) found that the dark matter halo 
can satisfy both the isotropy and inhomogeneity require- 
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merits of BATSE, however, a larger value of R corc than 23 
kpc is required in order to avoid having too few faint GRBs. 

In many halo models the nearest GRBs have a greater 
gradient, and thus, dipole moment than the entire distribu- 
tion so the anisotropy might be first detectable in a subset of 
near GRBs. This is not the case for the Uniform Spherical 
Halo since if we do not view far enough to see its edge, the 
model becomes isotropic and homogeneous. For the dark 
matter halo model, we find that for small data sets such that 
the total data set limit on R core is <R Q , examining the 
brightest \ or j of the GRBs gives a better limit. For larger 
data sets that already require a large value of R core , examin- 
ing a smaller data set gives a less strict limit because the 
fluctuations of the smaller sample size are more important 
than the increased dipole moment. In any case, the brightest 
sources, e.g., data sets 24, 25, 31, 32, 38, and 39, are consis- 
tent with isotropy. 

Smith (1994b) has proposed another geometric halo 
model, an exponential halo with source density oc exp ( — R/r). 
His example successful model has <cos 0} = 0.056, which 
deviates by 2.1 a from the value observed for 1005 GRBs. 
This particular exponential halo is a fair match for the data. 

5.6. High-Velocity Neutron Star Models 

The dark matter halo model with a large core radius can 
match the isotropy and M31 constraints, and also the inho- 
mogeneity constraint (Hakkila et al. 1994a). This model is 
primarily geometric because no candidate GRB source 
population has this distribution. This raises the question of 
whether the distribution of any known object can match the 
observations. Neutron stars are, of course, the most prom- 
ising candidates. 

Paczynski (1990) and Hartmann, Epstein, & Woosley 
(1990) calculated the orbits of neutron stars bom in the disk. 
They found that if the distribution was viewed deeply 
enough to detect a deficiency of faint sources there would be 
a significant quadrupole moment <sin 2 b - 3). New 
observations and analyses since then have demonstrated 
the existence of a large number of high-velocity neutron 
stars (see, e.g., Lyne & Lorimer 1994; Frail, Goss, & White- 
oak 1994). Lyne & Lorimer (1994) find the mean neutron 
star velocity to be 450 + 90 km s' 1 and the velocity dis- 
tribution to be very broad, with a root mean square width 
of 525 kms" 1 . The mean velocity is close to the local speed 
to escape the Galaxy, 550 ± 100 km s -1 (Fich & Tremaine 
1991). Many neutron stars are thus either unbound or only 
weakly bound, implying the existence of a much larger halo 
than formerly thought (Lyne & Lorimer 1994). However, 
these observations still require a substantial fraction of the 
neutron stars to have velocities below the Galactic escape 
speed. The low-velocity neutron stars generate a disk popu- 
lation in addition to the halo formed by the high-velocity 
neutron stars, so this model still has difficulty matching the 
observations. A variety of solutions have been suggested : (1) 
a hypothetical population of neutron stars that are born in 
the halo — this population either being more numerous than 
that born in the disk or more efficient at bursting (Brainerd 
1992; Eichler & Silk 1992; Hartmann 1992; Fabian & Pod- 
siadlowski 1993); (2) a delayed turn-on so neutron stars still 
near the disk are not seen (Li & Dermer 1992); (3) some 
mechanism such as a magnetic-field/recoil velocity corre- 
lation so that only high-velocity neutron stars produce 
gamma-ray bursts (Brainerd 1992; Li & Dermer 1992); and 


(4) an alignment between the recoil velocity of a neutron 
star and its magnetically beamed gamma-ray burst emis- 
sion (Li et al. 1994). 

Eichler & Silk (1992) propose that neutron stars are born 
in an extended halo via the mergers of Population III white 
dwarfs. They assume a dark matter halo distribution with 
Rcore = 40 kpc and estimate <cos 6 > to be 0.05, which devi- 
ates from the observed value by about 1.8 o. Hartmann 
(1992) calculated the orbits of neutron stars born with 400 
km s ’ recoil velocities in an R ~ 2 halo of radius 50 kpc. An 
observing depth of 30 kpc matches the intensity observ- 
ations but yields <sin 2 b - = -0.05, 5.0 o off from the 

current observations. Fabian & Podsiadlowski (1993) 
propose that gamma-ray burst sources are ejected from the 
Magellanic Clouds, producing an extended halo nearly 
homogeneous to 50 kpc. Their example model has a dipole 
moment to the Magellanic Clouds of 0.038, which deviates 
by 1.9 <t from observed moment (Table 8). The agreement of 
these models with the data ranges from satisfactory to unac- 
ceptable. 

Li & Dermer (1992) proposed that gamma-ray bursts 
originate after a time delay from neutron stars born in the 
disk with velocities of 1000 km s” 1 . The time delay allows 
many of the sources to leave the disk before they become 
active, thereby reducing the disk contribution. Their model 
has <cos =0.048, which is 1.7 0 off from the current 
observed value. 

The halo beaming model of Li et al. (1994) assumes that 
(1) only high-velocity neutron stars emit GRBs, (2) these 
objects are born with their magnetic axis aligned with their 
recoil velocity, and (3) the gamma-ray burst emission is 
beamed in a 20° half-angle beam along the magnetic axis. 
This beaming suppresses the observation of sources in the 
plane relative to those far above the plane. Li et al. (1994) 
determine the consequences of this model by following the 
orbits of neutron stars born in the disk with recoil velocities 
of 1000 km s *. This model still predicts a concentration in 
the plane, and, uniquely, a quadrupole concentration 
toward the Galactic center and anticenter; both concentra- 
tions should be most prominent for the nearest GRBs. The 
authors show, as a function of peak flux on the 1024 ms 
timescale, the model values of <sin 2 b - and <cos 2 6 
— 3) for sets of all GRBs brighter than the peak flux value. 
The brightest £ of the first 1005 GRBs deviate from the 
model values by 2.0 o and 3.1 o (Table 8). This particular 
form of the halo beaming model seems unlikely. 

The models of Podsiadlowski et al. (1995) use several 
methods to reduce the anisotropy of neutron stars ejected 
from the plane. They follow the orbits of neutron stars in a 
nonspherical halo potential, which randomizes the orbits, 
thereby producing greater isotropy. In order to avoid a disk 
signal, they assume that only high-velocity neutron stars 
burst and that the sources have a delayed turn-on and a 
turn-off time. The latter requirement prevents the build of a 
disk population formed by the capture of neutron stars that 
received their birth kick opposite to the Galactic rotation. 
Their example models (their Figs. 5a and 5b) are based 
upon: (I) standard candle sources observed to 340 kpc, (2) 
neutron stars born with recoil velocities between 600 and 
700 km s“‘, (3) a turn-on delay of 10 7 yr, and (4) turn-off 
times of 10 9 yr for model 5a and 10‘° yr for model 5b. The 
values of <cos 0> and <sin 2 b - }> of their model 5b are 2.0 
0 and 2.2 o off from those observed for 1005 GRBs, demon- 
strating that a source lifetime of 10 10 yr is probably too 
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long. The moments of the 10 9 yr lifetime model are in good 
agreement with the data, deviating by 1.4 a and 1.7 a. 

5.7. Two-Component Models 

It has frequently been suggested that GRBs might orig- 
inate from two populations in various combinations of disk, 
halo, and cosmological distances (see, e.g., De Jager & 
Schaefer 1992; Lingenfelter & Higdon 1992; Smith & Lamb 
1993; Higdon & Lingenfelter 1994; Katz 1994). Since suffi- 
ciently extended halo models agree with the observations 
we expect a small admixture of a disk component with 
either a halo or the isotropic distribution to also be satisfac- 
tory. As an example we consider a model with fraction / a „ is0 
of the locations drawn from the Watson distribution with 
k = — 1 and the remainder isotropic. The Watson distribu- 
tion with jc = — 1 has one-half of the probability within 
b 1/2 = 23? 1 of the plane and <sin 2 b — £> = —0.0796, an 
8.4 a detection for 1000 GRBs; an example simulation is 
shown in Figure lc. From Figure %g we see that BATSE 
would typically detect at the 99% confidence level a 30% 
contribution from this disk model in a sample of 1000 loca- 
tions. Of course, the isotropy constraints would allow a 
greater fraction of a less concentrated disk component. 
However, the contribution of an almost isotropic com- 
ponent is limited by its near homogeneity. 

Lingenfelter & Higdon (1992) consider a model with two 
Galactic components arising from a single source popu- 
lation. The low-luminosity emissions are seen to a nominal 
distance of 300 pc and the high-luminosity emissions to 100 
kpc. The source population is based upon the orbital calcu- 
lations of Paczynski (1990) for Galactic neutron stars born 
in the disk. The sources seen at the great distances contrib- 
ute most of the anisotropy — the low-luminosity sources are 
virtually isotropic and homogenous. This 1992 model now 
clearly disagrees with the data, differing from the observ- 
ations for 1005 GRBs at the 3.5 a level for <cos 6 > and at 
the 6.1 a level for <sin 2 b - £>. Their revised two- 
component model (Higdon & Lingenfelter 1994) is formed 
from two geometric models: an exponential disk and a dark 
matter halo model with a cutoff at 300 kpc. Since at most 
25% of the GRBs are from the disk, the dominant aniso- 
tropies of their models are those of the halo portions, the 
dipole moment <cos 0>. The comparisons are given in 
Table 8. The core radii of 7.5 and 15 kpc, which are 
“inferred from constraints on Galactic structure,” are in 
very poor agreement with the observations. Even the model 
with a halo component with a core radius of 30 kpc does 
not agree well with the observations. This last model has a 
disk component of only 8%. 

Smith & Lamb (1993) consider several other disk plus 
halo models. Their disk model is an exponential sampled to f 
of its scale height, resulting in a quadrupole moment <sin 2 b 
-£>=-0.0401, which is half the moment of our 
example two-component model (above. Figs, lc and 8 g). 
One of their models has 20% of the sources from the disk, 
while the remainder are from a dark matter halo with core 
radius 22.5 kpc observed to 135 kpc. Again, this is predomi- 
nantly a halo model, and it is currently 2.2 a away from 
BATSEs measurement of <cos 0>. Their other model has 
67% of the GRBs originating in the disk and the remainder 
from a Gaussian shell halo model. The halo model is a 
Gaussian distribution with <x = 38 kpc and its highest 
source density on a shell 25 kpc from the Galactic center. 
Since this model is dominated by its disk component, the 


crucial statistic is <sin 2 b - £> = -0.027, which is 2.6 o 
away from the value measured by BATSE for 1005 GRBs. 

In this paper we are emphasizing the isotropy character- 
istics of two-component models. It is also necessary for 
models to match the observed peak flux distribution, 
log N — log P . A particular requirement for two- 
component models to agree with BATSE's observations is 
that the slope of the logarithmic cumulative peak flux dis- 
tribution be concave down (Paczynski 1992; Smith & Lamb 
1993; Hakkila et al. 1994a; Harding 1994). 

5.8. Heliocentric Models 

No heliocentric pattern is evident (§ 4). The lack of a 
dipole moment toward the Sun is a weak constraint: a 
spherical shell at the 19 AU distance of Uranus would 
produce a dipole moment of <cos <£> = 0.035, which would 
be a 1.9 a detection for 1005 GRBs (adapting eq. [1] from 
galactocentric to heliocentric shells). The above limit 
applies to the “typical” gamma-ray burst; testing the 
plane-wave nature of the wavefront of bright and presum- 
ably close bursts yields lower limits of %40 AU (Conners et 
al. 1993). 

The models that have been suggested for a heliocentric 
origin involve comet-comet or comet-black hole collisions 
in the Oort cloud (Bickert & Greiner 1993; Katz 1993; 
White 1993; Luchkov 1994). One can, therefore, constrain 
these models based upon our knowledge of the cometary 
cloud surrounding the Sun: models of the distribution of 
the heliocentric cometary cloud are based upon a solar 
system origin, planetary and stellar perturbations, and the 
tidal force of the Galaxy (Weissman 1990; Clarke et al. 
1994). The Galactic tidal force has had negligible effect on 
the inner region of the cloud, which should, thus, have a 
strong ecliptic plane concentration reflecting its solar 
system origin (Clarke et al. 1994). The lack of an ecliptic 
plane concentration according to <sin 2 — £> is, thus, a 
useful constraint. In contrast, the distribution of the inter- 
mediate Oort cloud is dominated by Galactic tidal forces — 
Clarke etal. (1994) estimate <sin 2 b - £> = +0.12, which is 
quite inconsistent with the observed GRB locations. While 
the distant portion of the cloud is probably more isotropic, 
producing GRBs there but not in the inner portion is very 
implausible since the cometary density is lower and the 
comets are more distant (Clarke et al. 1994). 

A severe problem for Oort cloud models of GRBs is the 
homogeneity of the brightest GRBs (Maoz 1993b; Clarke et 
al. 1994; Horack et al. 1994), which has been well observed 
by instruments with greater exposure than BATSE (see, e.g., 
Fenimore et al. 1993). The physics of the models pose even 
greater problems (Clarke et al. 1994). 

6. DISCUSSION 

If gamma-ray bursts have a large-scale anisotropy, it is 
remarkably small, especially compared to any observed 
Galactic component. The most extended, observed com- 
ponent of the Galaxy is the halo, as traced by stars or 
globular clusters. The globular clusters are distributed in 
galactocentric radius R as R~ 35 (Zinn 1985; Gilmore, 
King, & van der Kruit 1990; Djorgovski & Meylan 1994). 
This contrasts with the more isotropic R~ 2 we used in the 
dark matter halo form (eq. [2]). Using a slightly different 
halo form than ours, Djorgovski & Meylan (1994) find a 
core radius of at most 2 kpc for the globular clusters, which 
contrasts with the 2 o lower limit for the core radius of the 



No. 1, 1996 


GRBs of 11 kpc. Zinn (1985) has shown that the globular 
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Galactic components. The observations currently require 
. a . 3 most a sma11 faction of the bursts originate from a 
sk component (§ 5.7). A common motivation for two 
population models is allow GRBs with cyclotron lines to 
originate from disk neutron stars in the solar neighborhood 
Models of photon emission from magnetic neutron stars 

ta7ce d to V CRB d T general,y believed that the dis- 

ance to GRBs was of order 100 pc. It seems unlikely that 

these models will work at halo distances that are >100 
times greater Two population models generally have a 
ratio of the distances of the populations of at least 100 
implying a ratio of luminosities of order 10 4 or more It 
seems implausible that phenomena with such greatly Se - 
en it luminosities would not be observationaJly separable 
nto two classes. The observed bimodality in GRB dura- 
tions (Mazets et al. 1981; Hurley 1991; Klebesadel 1992- 
Dezalay et al. 1992; Kouveliotou et al. 1993b) is the best 
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distnbut.onswith the observed log N - log P distribution 
The simultaneous observations by BATSE of the isot- 
ropy and the strong inhomogeneity of GRBs disprove the 
previous paradigm of a Galactic disk origin. No known 
Gaelic component, including the Dark Matter halo is 
sufficiently isotropic to match the observational constraints 
Consequently, current models either postulate a new Galac 

Uo3r‘, 0 H eX,reme ' ner ® firemen, s ?o, 

SKi 'T? GRBs are GaIac,ic - ,he lack of 

excess ORBs from the Galactic center and M31 constrain 

k 1? d,slan “ of the typical GRB observed by BATSE lo 

Srfhe Z and 8 5 5k conlras t .O the difficulty 

that the isotropy observations pose for Galactic models 

otropy is a natural consequence of an origin at cosmo- 

Ws bufdJT 8 ' T hC ° bserv . ed lar 8 e ' sca le isotropy strongly 
lavors bu does not require that gamma-ray bursts originate 

at cosmological distances. The continued observation of 

gamma-ray bursts by BATSE promises to either detect an 

models° Py ° f ‘° fUrther ‘ ighten the constrai nts on Galactic 
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